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ABSTRACT 

This document describes a theoretical research program on the undercooling 
and solidification of materials under variable external field conditions 
performed by The Boeing Aerospace Company for the National Aeronautics 
and Space Administration. Work was performed under the direction of 
Dr. Mary Helen Johnston of Marshall Space Flight Center, the Contracting 
Officer's Representative. 

A catalog of theories and models of nucleation of solid phases in the 
melt is provided in the document, as is discussion of the relation of 
undercooling to intermolecular potentials, the dependence of grov/th rate 
on undercooling, the influence of undercooling on liquid^solid interface 
stability and sol id structure, the direct effects of external fields 
on melts, the relation of solid physical properties to structure and 
the role of nucleants in solidifica‘:ion. Results of theoretical analysis 
of several experiments related to the Space Processing Applications 
Program are given in the document, and recommendations for future 
experiments and further theoretical development along v/ith procedures 
for correlation of theory and experiment are specified for solidification 
in space research. 
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PART 1 - SUMMARY 

1.1 STUDY OBJECTIVES 

It 1s commonly accepted in science and engineering that a phenomenon 
is not really understood until a mathematical theory of the phenomenon 
is developed from which quantitative predictions can be made which 
agree with experimental studies of the phenomenon. Thus the over-all 
purpose of the study "Undercooling of Materials During Solidification 
in Space" is to theoretically determine how undercooling, gravity and 
nucleants affect the solidification of materials from the melt and 
influence the properties of the resultant solid. Specific objectives 
which complement this purpose include: 

1. Gain an understanding of the theoretical relationship between 
undercooling and both homogeneous and heterogeneous nucleation. 

2. Determine the influence of external fields (particularly gravity) 
on the nucleation process. 

3. Gain an understanding of how solidification history affects the 
physical properties of a material. 

4. Theoretically analyze some solidification experiments performed 
at MSEC and as part of the Space Processing Applications Rocket 
(SPAR) program. 

5. Specify materials to be used in the comparison of theory to 
experiments in one-g and zero-g. 

The twelve-month research program described by this report constitutes 
a major step toward achieving these ambitious objectives. 
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1.2 APPROACH 

The basic approach to this study vms to obtain knowledge of the theory 
of nucleation from the melt (nucleation of solid or liquid phases in a 
vapor phase was not considered) , and adapt the theory to Space Pro- 
cessing program needs and experiments in which the acceleration 
due to gravity, g, is a variable. In cases v/here Space Processing 
experimental situations did not correspond to existing theory, original 
theoretical relationships were developed to carry out the required 
analyses. 

The above approach was organized around four tasks: 

1. evaluation of Nucleation Theories 

2. Influence of Undercooling and Nucleation on Properties of Solids 

3. Investigation of Field Effects on Undercooling and Nucleation 

4. Correlation of Theoretical and Experimental Results 

The theories and models identified in Task 1 were used in Tasks 2, 3 
and 4 to explore previously unanalyzed experimental results, and to 
begin development of a theoretical model to predict solid material 
properties from the melt properties and the solidification environment. 
Environmental variables considered include gravity level, nucleants 
(impurities, particles of high temperature phases and container 
surfaces) and temperature (undercooling). 

Theories and models of nucleation from the melt were obtained from the 
literature, and included both classical and modern homogeneous nucleation 
theories, heterogeneous nucleation theories and models and energetic 
aspects of nucleation in tv;o-component melts, [‘•lost of cl assicaT nucleation 
theory is based on the work of David Turnbull. His experimental and 
theoretical results form the basis for much of this study (see References). 
These theories were used to derive equations which show how undercooling 
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1*2 (Continued) 

influences liquid-solid interface stability, growth rate, nucleation 
rate and structure of material solidified from the melt. An approach 
to relating solid physical properties to structure was developed for 
use in the comprehensive model for solid property prediction. Other 
theoretical research included: an explanation of why some pure 

materials sustain higher undercooling (before homogeneous nucleation 
initiates) and others ; how nuclei segregate in a melt under the 
influence of acceleration fields; the direct effects of external fields 
on solidification; and some mathematical analysis of experiments 
performed at MSFC, in SPAR flights, and experiments performed by 
Kattamis, et.al., v-zhich are applicable to Space Processing goals. 

1.3 RESULTS OF THE STUDY ' 

Hucleation Theory 

A catalog of theories of nucleation from the melt has been created 
(see Section 2.1), containing the theories of classical and modern 
homogeneous nucleation, heterogeneous nucleation and nucleation 
in two-component melts. For both homogeneous and heterogeneous 
nucleation, models are given v/hich describe the formation of nuclei 
of different shapes. The most common case, heterogeneous nucleation, 
is most easily described in terms of a spherical cap shaped nucleus 
and the angle the tangent to this cap nucleus makes with the nucTeant 
surface, called the wetting angle, 0 . This angle is related to the 
interfacial energies betv/een nucleus and melt, nucleus and nucleant 
and nucleant and melt. Other parameters related to these same inter- 
facial energies are derived for different nucleus shapes, but since the 
exact shape of a nucleus is seldom known, e provides a convenient 
parameter for characterizing this aspect of nucleation. Nucleants 
considered included particles of high temperature phases suspended in 
the melt as well as other impurities which might be present, and 
melt container surfaces. A major question in nucleation theory, which 
remains unanswered as yet, is the exact mechanism by which a particular 
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1.3 (Continued! 

nucleant causes a nucleus to form. Research reported in the litera- 
ture centers on disregistry between the lattice structure of the 
nucleant and that of the nucleus, but there is no agreement among 
authors as towhat the mechanism involved is or how to describe it 
theoretically. 

f^aximum Undercooling 

The question of why tv/o different pure materials, with no nucleants 
present, will homogeneously nucleate at different undercoolings has 
been resolved by showing that the maximum undercooling a pure melt can 
attain 1s proportional to the ratio where ct is the liquid- 

solid interfacial energy and S is the entropy of fusion for the 
material (see Section 2.2.1). is related to the intermolecular potentials 
of both the liquid and solid phases by statistical mechanics, cr can be 
related to liquid or solid phase intermolecular potentials by several 1 

semi -empirical theories, most of v/hich involve structural models of 
the liquid-solid interface. Thus the maximum undercooling attainable 
in a pure melt is a (complicated) function of the intermolecular 
potentials of the material, and hence is a unique characteristic of 
that material . 

Environmeiital Effects on Nuclei 

k 

Experiments by Dr. Johnston at MSEC have shown a pronounced solute segre- 
gation effect in lead-tin eutectic material as a function of the environmental 
variables gravity level (acceleration) and undercooling (temperature). 

Although an increase fn segregation with gravity level v/as expected, a segre- 
gation increase with increasing undercooling yras not. To analyze this situation, 
it was assumed that nuclei settle in the melt due to gravity, and segregate 
further due to an unexplained undercooling effect which may or may not be 
coupled to the gravity effect. The gravity effect was analyzed by solving 
the equation of motion for a nucleus moving in the melt and substituting ^ 
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1.3 (Continued) 

the result into Pick's law for concentration distribution in the melt 
(see Section 2.2.4). The resulting differential equation can be solved, 
yielding a transcendental equation for solute concentration as a function 
of position in the melt, gravity and undercooling. This equation is in 
good agreement v/ith experimental (one-g) concentration distribution 
data reported in the literature. It also shows qualitative agreement 
with the MSFC experiment for concentration versus gravity level, if 
allov/ance is made for the unexplained undercooling effect. There is 
no agreement between theory and experiment for concentration as a 
function of undercooling, the experimental data being approximately 
linear and the theoretical equation being strongly exponential. 

Several theoretical hypotheses for this undercooling effect were 
tested, but none showed agreement with the observed experimental data. 
Further research is needed to identify the mechanism of this anomalous 
undercooling effect. A paper detailing the above v^ork was presented at 
the Winter Meeting of the American Physical Society, held at California 
Institute of Technology (Bull .Am.Phys.Soc. , 1511 , 1975). 

Solidification Model Development 

As part of the over-all goal of this study to ascertain how solidification 
environment affects resultant solid properties, an attempt was made to 
develop a mathematical model which would predict grain size and dendrite 
volume, morphology and spacing as a function of nucleant w'etting angle, 
e, and the undercooling, AT (see Section 2.2.7). This model, in its 
present form, solves 8 simultaneous equations; those for nucleus volume, 
dendrite growth volume, dendrite growth rate, number of nuclei created, 
micleation rate, nucleus critical radius, maximum solute diffusion time 
and maxi mim heat flux time. A Wang 700A computer program was written 
for calculations with this set of equations. The model has been compared 
to experimental studies by Kattamis and Flemings on iron-nickel alloy, 
and shoves rough agreement with the. experimental data over limited ranges 
of undercooling. But for the most part, the theory does not yet reflect 
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1.3 (Continued) 

the experimental findings^ largely because of simplifying assumptions 
which lead to rather unrealistic equations. A great deal of additional 
time and effort will be required to derive rigorous expressions for such 
parameters as total solidification time, dendrite growth volume and 
solute diffusion time. Model development thus far does, however^ provide 
a firm basis for further research, and also provides grounds for the 
hope that the solidification model can be refined to give reasonably 
accurate predictions of grain size and dendrite volume and morphology 
data. 


Comprehensive Property Prediction Model 

Work v/ith the solidification model has lead to the concept of the 
comprehensive property prediction modal (CPPM) (see Section 2.4.1), 
which is to be composed of several elements (comprehensive refers to 
its multi-model composition), each of which models a different part 
of the solidification process. Output of the CPPM would be solid 
material properties calculated from input consisting of melt property 
(constant) data and external environmental variables, such as gravity 
level, nucleant wetting angle and cooling rates (v;hich lead in the heat 
transfer model to undercooling values). Once the solid properties have 
been output, they can be compared to material property requirements 
specified by the potential user of the material. If calculated values 
are too far from the required values of material properties , new 
melt properties (signifying different raelt compohents or solute 
concentrations) along v/ith altered environmental variable ranges, if 
needed, can be input again, and this iteration repeated until a match 
between calculated solid properties and user requirements is obtained. 
In addition to its usefulness as a molecular-level materials design 
tool, the CPPM could be used by the Space Processing Applications 
(SPA) program as a guide for the definition and design of SPA 
experiments and as an aid in selecting the most promising Tf:ater1al 
combinations (from among the 100 million or so two-to-five component 
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1.3 (Conti’iued) 

alloys possible in microgravity) for flight experiments. Therefore, 
development of a CPPM should receive emphasis v/ithin the SPA Program. 

Relation of Physical Properties to Structure 

One element of a potential CPPM which received attention during this 
study was the theoretical relationship between material structure 
and physical properties (see Section 2.2.6). Most solid physical 
properties derive from the periodic nature of the crystal lattice and 
its symmetry with respect to spatial operations and the electron and 
phonon energy spectra characteristic of the material. However, a few 
physical properties are influenced by non-periodic structures in the 
solid, such as grain boundaries, impurities or other sites of electron 
and phonon scattering- In particular, three physical properties of solids, 
the thermal conductivity, the electrical conductivity and the thermo- 
electric power have been found to depend on the phonon mean free path 
in the solid. The phonon mean free path is itself related approximately 
to grain size in the solid. The dependence of these properties on 
grain size is usually masked by the fact that thermal conductivity is 
proportional to heat capacity as welT as phonon mean free path, and 
heat canacity varies much more rapidly vnth temperature than mean free 
path varies with grain size. Thus thermal conductivity, electrical con- 
ductivity and thermoelectric power are influenced by grain size in the 
solid, though to a lesser extent than by temperature. 

SPAR Experiments 

Two SPAR experiments to investigate the solidification of the metal- 
model material , anmioni urn chloride solution (NH^Cl-H^O), were flown 
during the study period (see Section 2.3.1). Certain aspects of these 
experiments were analyzed theoretically: the undercooling of the 
first sample corresponding to measured growth rates, and the failure: 
of the second sample to solidify, even v/ith substantial undercooling. 

OKTGMAi; PAGE IS 

OF POOR QUALITY 
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Very fev/ of the melt property parameters for l'!H^C 1 -H 20 are known, so 
to do calculations for this material, nominal metal values (since it is 
a metal -model ) were assumed for entropy of fusion, S , and diffusion 
coefficient, and the liquid-solid interfacial energy, a, v/as estimated 
with a relationship between a and due to Turnbull. These properties 
were then used to compute values of growth rate versus undercooling 
for the Flight 1 experiment and nucleant wetting angle versus undercooling 
for the Flight 2 experiment. From these values it was learned that the 
first sample did indeed undsrcool by an amount between 1®K and 4'^K, and 
that superheating of the melt prior to launch of the second sample 
dissolved or “denucleated" all nucleants characterized by wetting angles 
up to 62®. Denucleation of the Flight 2 experiment had already been 
advanced as the reason for its failure to solidify. The theoretical 
analysis provided a quantitative description of the denucleation » 

F11 ght Experiment Materi a1 s Sel ecti on 

One task of this study was to select two material systems for micro- 
gravity experimentation for the ourpose of comparison of theory and 
experiment. A major problem in comparing theory to experiment is the 
unavailability of established values for the liquid state properties 
of materials - especially alloys. There is virtually no data on the 
melt properties of multi -component materials, and very little data on 
melts of elements. Even though there is no good method of calculating 
alloy melt properties even when the elemental component property 
values are known, it is felt that theory-experiment comparison could 
best be carried out with materials systems for v/hose component elements 
pertinent melt properties are known. In addition to the availability 
of melt property data, the following criteria influenced the selection 
of material systems for flight experimentation; 

9 Abundance of previously published ground-based research on the 
material system which may prove pertinent to Space Processing 
■objectives...; 
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1.3 (ConUnued) 

0 The practical (econDHiic) importance of gaining more knowledge about 
the material system, 

9 The e/.tent to which the material system properties improve the ability 
to study a particular phenomenon, and 


0 The general simplicity of the material system with respect to 
experiment design and theoretical analyses. 

Based on these criteria, the material systems selected are iron- 
nickel and lead-tin. The melt properties for these elements are given 
in Table VII. . 

1.4 CONCLUSIONS AND RECOMMENDATIONS 


The major conclusions of this study are as follows: 


1. Classical heterogeneous nucleation theory is at present the 
simplest method for performing calculations of solidification phenomena. 
The validity of some assumptions in the theory have yet to be established, 
and its accuracy is impaired by a lack of values for the melt properties 
which comprise the input to the theory. But these problems also exist 

in the modern theories of nucleation, which are conceptually and 
computationally more complicated. 

2. Gravitational settling (segregation) of dense phase nuclei in a 
melt is described by the equation 


B,C + ^1 


XI) 

B ^ 


^ yC .V'/B, 


0 



where C is solute, concentration, g is acceleration due to gravity and 
z is axial position in the melt. /is a parameter which depends on 
undercooling. The other factors are combinations of melt properties 
which are constant for a given material . 
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1.4 (Continued) 

3. The segregation of lead and tin dendrites observed in Pb-Sn eutectic 
material ei<perirnents in the MSFC centrifuge, can be only partly explained 
by the settling theory discussed above. The linear dependence of 
concentration on undercooling, at fixed 2 , is as yet unexplained. 

4. The maximuin undercooling attainable in pure melts is determined 

by the intermolecular potentials in the liquid and solid phases of the 
material . 

5. The amount of undercooling achieved in the December 1975 SPAR flight 

experiment with ammonium chloride solution was predicted by theory to be 
between 1°K and Classical theory also predicts that the May 1976 SPAR 

ammonium chloride flight experiment ms denucleated to the extent that 

all nucleants characterized by wetting angles of up to 62® were dissolved 
into the melt. 

6. There 1$ strong theoretical and experimental evidence that grov/th 
rate of a solid into its melt is a function of the undercooling squared. 

7. Physical properties of solids can be related to grain size through 
phonon scattering theory, but mathematically, the dependence of these 
properties on grain size is weak compared to their dependence on 
temperature. 

8. It should be possible to develop a solidification model which would 
predict grain size and dendrite morphology as a function of undercooling 
and nucleant characteristics. The major obstacle appears to be the 
lack of melt property and nucleant characteristic data for input to 

the model. It may also be possible to develop a comprehensive material 
property prediction model which v;ould facilitate the design of materials 
to meet specific, predetermined requirements by calculating the material 
composition and solidification environment (gravity level , temperature 
profile, container type or absence) needed to produce the desired 
material . 


1-10 
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1.4 (Continued) 

9. Two material systems v/h1ch are well suited to the comparison of 
theoretical predictions to experimental reality and to the objectives 
of the Space Processing program are the systems iron-nickel and lead- 
ti n . 

Based on these conclusions, it is recommended that there be further 
study of: 

a. Computational methods and their accuracy within nucleation theories, 

b. The relation of undercooling to solid structure, including the solute 
segregation effects observed at M5FC, 

c. The relation of solid structure to all properties of the solid, 

d. The relation of liquid, solid and nucleant structure to the 
nucleant parameters { 0 , Aa) and to undercooling. 

Additional recommendations are: 

e. That experimental measurements be made in microgravity of undercooling, 
nucleation rate, grov/th rate, grain size and morphology in the Pb-Sn 
and Fe-Ni systems, 

f. That a program be initiated for the measurement and/or calculation 
of liquid state (melt) properties of engineering materials, such as 
the Fe-Ni alloy system, 

g. That development of a comprehensive material property prediction 
model be undertaken. Such a model would be used as a guide for 
Space Pv’ocessing experiment definition and design, as a guide for 
selection of materials for Space Processing experiments which 
show promise of commercial value, and as a molecular-level 
materials design tool. 
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PART 2 - DETAILS OF STUDY RESULTS 


2.1 THEORIES AND MODELS OF KUCLEATION 

The solidification process is possibly the least well understood 
phenomena of those which are common to everyday life. This is demonstrated 
by the plethora of theories and models which have been developed to 
attempt a quantitative description of the process, and by the failure 
of most of these models to explain any large body of experimental 
observation. These models as a group do, hovjever, provide significant 
insight into solidification phenomena, and a specific model properly 
applied over a limited range can yield quantitative predictions which 
agree quite well with the experimental data it seeks to explain. Also, 
the theories and models about to be discussed are the best available 
mathematical tools for making predictions about nucleation. This 
study is concerned with solidification from the malt, so theories 
of nucleation of a condensed phase from vapor will not be considered 
explicitly, although much of the mathematical formalism is common to 
both nucleation events. The basic goal of most nucleation theories 
is to calculate a nucleation rate of the form 


AGr 


I - 




kT 



( 1 ) 


where 

is a frequency term (nuclel/time/volume) 

aG„ is an activation energy term for diffusion across a liquid-solid 
D 

interface 

aG* is the free energy required for formation of a critical nucleus. 

The different models are required to calculate the terms AGp and 
aG* for different physical situations. 
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2.1,1 Homogeneous Hucleation - Classical Theory 


Classical nucleation theory has been derived from absolute rate theory 
by Becker and Coring [1], Becker [2}, Turnbull and Fisher [3], Turnbull [4], 
Hollomon and Turnbull [5] among others [6]. Homogeneous nucleation refers 
to the spontaneous formation of solid nuclei in the melt of a pure 
material. A spherically shaped nucleus will be assumed to demonstrate 
the classical approach. 


The free energy of formation for a spherical nucleus is 

aG = " Trr^ AGy + 4irr^a (2) 

where AG = S AT is the bulk fv’ee energy with S being the entropy of 

V L# • 

solidification and AT the undercooling (T^^ - T). a is the surface energy 
at the nucleus-melt interface. Maximizing ag with respect to nucleus 
radiusa r, yields the expression for the critical radius 


r* = ||~ . (3) 

Substitution of (3) in (2) gives the critical free energy 

3 

AG* = AG(r*) = y . (4) 

3S^^aT^ 

If there is -reason to assume nuclei of different shape, the procedure 
remains the same, but models of different nuclei shape yield: 


Cylinder: r* = 


2cr . ^ _ 4g . p* _ „ cr 

p a h ,p a Au — 0(1 o 

(AGJ^ 


4cr 4a. 

Parallelepiped: a* ~ ~rf~ ^ Tr~^ ^ 


* = 


4a 

c 


ag* 


32 


a b c 
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2.1.1 (Continued) 


Several attempts have been made to determine and aGj^. Some of the 
more useful results are; 

0 Turnbull and F1sher[3]: - n‘-^C ~ 

where n* is the number of atoms in the surface of a critical 
nucleus, 

C is the number of atoms per volume in the melt, 
h is Planck's constant, 

and a is a parameter dependent on nucleus shape. 

0 Free volume model [7]: re 

where r^ is the radius of an atom of the subject material. 

In practical calculations of nucleation rate, however, I is so sensitiv 
to exp(- AG*/kT), that the remaining terms may be approximated by a 
constant = I' exp{” AG^/kT). It has also been shown [4] that for 
several material systems, I^ may be taken to be 

T - in33 „ 75.985 

1 “ 1 U “ e 

• ■0 

with very little error in the final value of I [6], So this is 

the value used for I in this document. 

0 . . . 

2.1.2 Homogeneous Nucleation -Modern Theory 

Modern homogensous nucleation theory attempts to calculate the free 
energy changes associated with nucleus formation (whatever the mode of 
formation) from statistical mechanics. This alv/ays involves the 
computation of a partition function, Q, since the free energy of a 
material system in a condensed state is derived from the partition 
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2.1.2 (Continued) 
function by 

G = F = -kT In Q (6) 

where F 1s the Helmholtz free energy. 


Replacement Partition Function 

An important portion of modern nucleation theory deals with the calculation 
of a "replacement" partition function, \gp- 

The partition function for a niicleus containing n molecules in a melt 
containing C(n) nuclei per unit volume is [8] 


Q = [- 


ClnT! 


■] 


1 

c(n) Xq, 


exp [‘ 


'rep 


aA(n) 

kT 


nu~ ] 

|.T-i 


(7) 


where X ~ h"^ {2irnmkT)^^^ is the translational partition function, 


Qj, “ h~^(2kTI)^'^^ is the rotational partition function, I 

being the moment of inertia, nmr^, of the nucleus, a is the 
liquid-solid interfacial energy, A(n) is the surface area of 
the nucleus, and u is the chemical potential of the bulk solid. 


The replacement partition function arises from the following consideration. 
When a nucleus forms from a melt, the atoms in the nucleus give up degrees 
of freedom v/hich are replaced by six translational and rotational degrees 
of freedom for the nucleus as an entity in the melt, with the relative 
positions of the atoms in the nucleus remaining fixed. That is, q^^p 
yields correction terms for the thermodynamic functions which are 
necessary to extend the validity of these macroscopic functions to non- 
macroscopic systems suclv as nuclei . Hence the name "replacement" 
partition function. Another v/ay to view this process is that as the 
atoms arrange themselves in the ordered structure of the nucleus, the 


2-4 



D256-1 0202 


2.1.2 (Continued) 

entropy of the system is lowered^ giving rise to an increase in free 
energy » calculated from 

rep 'rep 

Crystal Subdivision 


One approach followed in modern nucleation theory is to calculate the 
change in free energy of a nucleus due to a change in surface area from 
partition functions for the system before and after the change. Part of 
this approach considers extension of a crystal surface^ while another 
part considers the energy change when a crystal is subdivided. In the 
latter part, the free energy per crystallite (v/hich divided itself 
from the parent crystal of n atoms) is 


3x-6 1 3n-6 


( 9 ) 


where there are m nuclei containing x molecules for a total of n = mx 
molecules in the solid, and where 


e"A 
X x 


increase in potential energy of a nucleus with surface area 

V - 


= translational and rotational partition functions, 

J.I- 4*h 

“ free energies of the and j ^ oscillation modes, 

th th 

,tij == zero point free energies of the i and J modes.. 

The i subscripts refer to the bulk, undivided crystal , while the j 
subscripts refer to internal vibrations of each crystallite. In terms 
of the classical surface free energy, o = e - Ts_, this is: 

^ CO CO w 


MGIKAnPAGEIS 
m POOR QUALOT 2-5 
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2.1.2 (Continued) 


AG^ = - kT q^) - kT 1 * kT In (10) 

V'/here 


kT In K ' 




3x-6 

i: 

j-1 


(♦j+nj) 


1 

ni 


3n-6 

s (f, -+£,•) 

i=l 


1 


(11) 


Evaluation of kT In !C is a central problem in modern nucleation theory. 

A 

This has not yet been done rigorouslyj but some prescriptions have been 
given for approximate evaluation. [9] 

Q Frenkel Prescription: kT Zn ~ ~ 3 ^" 


Note: the subscript ® refers to a flat surface. 


0 Kuhrt Prescription: kT Zn 


Note: Ujjco “ crystal chemical potential. 

0 Lothe-Pound Prescription: kT Zn K* = {o -a ) A -i- 2 

X ^ 


Ts 


X »V X 


2.1.3 Heterogeneous Nucleation Theory and Models 


Heterogeneous nucleation refers to the formation of a nucleus of some shape 
on a foreign substrates whether this substrate is a container surface ^ 
the surface of a contaminant particle or the surface of a particle of 
a high temperature phase of the melt material itself. Equation (1) 
holds for heterogeneous nucleation as well as for homogeneous nucleation 
if the appropriate expressions for AG^ and aG* are used. Equation (2) 
is generalized as follows 

AG = '’n(p^-y^) •!' a A (12) 

where n = number of moles in a nucleus 

= chemical potential of liquid and solid respectively 
a = interfacial energy 

A = area of liquid-solid interface (any shape). 

2-6 ...... 
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2.1.3 (Continued) 


The difference in potential, yp-i-i-# is called the driving force for 

S 

solidification, and the volume potential is defined by 







(T3) 


where V is molar volume. The derivation of the classical example of a . 
spherical cap nucleus v/i 11 be given, then the results of critical 
parameters for other geometries will be quoted. 


Spherical Cap Nucleus 

Figure 1 shov/s a schematic of a 
spherical cap nucleus, n, on some 
kind of substrate, s. The cap's 
radius is r, and the interfacial 
energies between liquid and sub- 
strate, £s, liquid and nucleus, 
nC, and nucleus and substrate, sn. 



are as shown. {Actually, interfacial 
energy, a, is a scalar quantity, 
and the vectors shown are of 


FIGURE 1 


SUBSTRATE. 

(IMPURITY): 

SPHERICAL CAP MUCLEUS 


interfacial tensions, which are 
vector quantities; but in most 
cases the magnitude of an inter- 
facial tension is numerically equal to the interfacial energy, particularly 
if there is little or no temperature dependence of these quantities). 

In this case, equation (12) becomes 


AC* = - nV r 


^ TOGS- 


PAGBlS 

QOAlSEfi 


(H) 


Z-7 
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2.1.3 (Continued) 

If R is the radius of the cap nucleus-substrate interface, then from 
Figure 1 


ttR^ = Trr^sin^o - 7rr^ (1-cos^€) 

since R = i" cos (|- - e) = sin e. Standard mensuration formulae may 
be used to find the surface area and volume of the nucleus 

2 

A-j = 2Trr (r-r cos e) = 2-Err ( 1-cos e) 

V = (r-r cos [3r-(r-r cos o)] = r’^ (2+cos e) (1-cos e)^ 

or 


V = Itt r^ f (e) 


1 2 

where f(e) = ^ (2 + cos 0 ) (1 - cos 0 ) . Thus equation (14) becomes 


aG = - ■|7rr^f(0)n ^\i^ -t- 2irr^(l-cos e)a^jj,+7rr'^(1-cos^o)(a2^“Cf2p) . 


2/1 2 
(1-cos 

This equation can be simplified by noting from Figure 1 that 


^sn " 


- cos 6 


Thus, after performing the algebraic manipulation, one obtains 
aG = - -|nr^f(0)nAUy + 4irr^a^j^f(9) . 

Haximizing aG with respect to r in the standard fashion 


(15) 


06 ) 


07 ) 
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2.1.3 (Conti luted) 


Thus the critical (sphere) radius per mole of the nucleus 1s 


* „ na lU 

Ay S aT 

C 

v/hlch Is the same as for the homogeneous case, 
one obtains for the critical free energy 


(19) 

Substituting (19) in (17)> 


ISti 

aG* = ^ f{8) (20) 

3Ah/ 

Which differs from the homogeneous case only by the "heterogeneous 
nucleation function," f(e) where e is called the "wetting angle". 

Models of Other Nucleus Geometries 

To discuss other geometries it is convenient to define the "heterogeneous 
nucleation parameter," Aa, by 

= i^ni. - “ds + “sn ' 


Related to the spherical cap model, this is just 


2act = (1 - cos e) 


(21b) 


However, in other nucleus geometries where 0 has no meaning, the critical 
parameters may be written in terms of the heterogeneous nucleation 
parameter as is done in these two examples [10]. 


Cylinder: r 


2c 


Ay 


ns, ^ 

— , h" " T — 


AG* - Sir 


2 ^ 
U»v)‘ 


n. 4a. 4o 

Parallelepiped: a^^ = 7 —^ > 77 "’ " 7i — => ~ — I" ’ 

% : : (Ay^)^ 
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2J.3 (Continued) 

Kossel Crystal 

One of the simplest heterogeneous nucleation models is the Kossel crystals 
shown in Figure 2. The basic assumption in this model is that the energy 
of attachment of a molecule to the nucleus 
is given in terms of the nearest-neighbor 
interaction energy, It is also assumed 
that, since the nucleus is orthorhombic 
and the molecules fit into a rectangular 
lattice, the interaction energy is related 
to the interfacial free energies by 


A = 2b c a ; d). = 2a c a, ; A = 2a b a 
^a 0 0 a’ ^b o o b’ o o ( 



where the ‘'a'‘face is that face orthogonal to the "a" direction in the 
crystal and bounded by sides b and c. The a^, b^ and c^ are the dimensions 
of a single molecule- in the lattice. Since the free energy of formation 
of a nucleus is in general 


aG = z - An^V 
i 

for the rectangular nucleus we have 


( 22 ) 


AG = 2bca, 2aca, + 2aba^ - abc A p., 

3 D C V 

The critical values of the nucleus dimensions and the free energy are 
found by maximizing aG with respect to a, b and c, and are 


4cj 


Ay. 


, b'^ = 


4aj 

Ay, 


c^ = 


Ay, 


and 


iG* = 32 


2 

(iPy) 


(23) 


(24) 


(25) 
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2.1.3 (Continued) 

Nucleus-Substrate Lattice Disregistr.y 

In earlier work [11] on undercooling of multiphase alloy systems* it v/as 
suggested that tiny particles of high temperature phases extant in 
the melt could act as heterogeneous nucleation sites for the lower 
temperature phases and prevent their undercooling. The theoretical 
approach to the problem of determining v/hether particles of a high 
temperature phase act as ’’nucleation catalysts” would be to ask first 
v/hether high temperature phase particles are more effective as 
heterogeneous nucleants than other Impurities in the melt. If the 
answer to this question is yes, high temperature phase particles are 
more effective nucleants, the next question is, why? 

The original work (pertinent to this question) by Turnbull and Vonnegut 
[12] defined a disregis try parameter, 6 by 

6 = AS. (26) 

0 

where a^ is the dimension of the unit cell of the solid material being 
nucleated and Aa is the amount by which the dimension of the nucleant 
(suostrate) unit cell exceeds a^- Turnbiill and Vonnegut concluded from 
the data available at that time that the smaller the disregistry, 6, the 
more effective as a nucleant was the particle or substrate catalyzing 
nucleation. They even shov/ed that the free energy of solidification 
required to produce a perceptible nucleation rate v;as proportional to 
6^. However, subsequent experimental v/Ork by Sundqui st and Mondol fo [13], 
essentially disproves Turnbull and Vonnegut' s theory showing that: 

(a) substrates or impurities with high disregistries are as good nucleating 
agents as substrates with low disregistries, (b) symmetry of the lattice 
plane orientation (across the mating plane) is at least as- important as 
disregistry in catalyzing nucleation, and (c) neither the orientation 
relationship nor the disregistry appear to have any effect on the under- 
cooling (and hence driving force, i.e., Gibbs free energy) required for 
nucleation. Yet work by Glicksman and Childs [14] at almost the same time 

Z-n y, 
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2.1.3 (Continued) 

found a qualitative V’elation between disregistry and undercooling;, more 
or less in agreement with Turnbull and Vonnegut. Hence the question of 
the influence of substrate (or impurity) surface structure on the value of 
undercooling at which nucleation occurs remains unanswered, and it must 
be answered before progress can be expected on determining the influence 
of high temperature phase particles on the melt. In fact, experimental 
studies of certain alloy systems show that neither oxides nor carbides 
were as effective as nucleants for lower temperature phases as might 
be expected, nor v/ere they as effective as particles of unknown contaminants 
present In the melts of the lower temperature phases [6], This v/ould 
indicate that concern over the effectiveness of high temperature phase 
particles as nucleants is no more justified than for any other contaminant 
in the melt. But there is no doubt that the lattice structures at the 
nucleus-substrate interface are important for determining the effectiveness 
of the substrate as a nucleation catalyst, even if the structure effect 
should manifest itself only in the interfacial energies shov/n in Figure 1. 
This is one area of Materials Science where a major research effort is 
urgently needed. It is possible that the theoretical portion of this 
research should proceed along the lines to be suggested in Section 2.2.1 . 

2.1.4 Hucleation in Two Component Melts 

Nucleation in two-component melts can best be understood with the help 
of the free energy-composition diagram shown in Figure 3. The parameters 
to be calculated are the criticaT radius of the nucleus, r*, and the 
critical free energy for nucleus stability, aG*. The change in bulk free 
energy upon formation of a nucleus can be expressed as the sum of the 
free energies after formation minus the melt free energy before 
formation [15] 

= % of atoms in nucleus X G of solid 

■t* % of atoms in depleted melt X G of depleted melt 
- G of original melt 
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2.1,4 (Continued) 


or, as can be seen from Figure 3, for a nucleus of composition 
which solidifies from an original melt of composition X leaving a 
depleted melt of composition X-j, 



0 X’ X Xj 1.0 

Xb 

FIGURE 3 GIBBS FREE ENERGY OF A TWO-COMPONENT MATERIAL 


Here + Ng “ v/here number of atoms in the nucleus, N-j is 
the number of atoms in the depleted melt and N is the total number of 
atoms in the binary system. Thus (27) can be rewritten 




Gj,) + 






The following compositional relationship holds: 


# nucleus atoms [nucleus composition - pure melt composition] - 

depleted melt atoms [pure melt composition - depleted melt 
composition] 
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2.1.4 (Continued) 
or 


Ng (Xg " X) = (X - X^) 
which is 


N. X - X. ■ • 

(C " X-, - X • 

* N, li. 

Substituting (29) in (28), by noting that 


yi el ds 


(29) 


X - Xp 

" IT X^' ^ X 




^X. " ®x 

Expressing - p z '~% — ^ differential instead of a difference, (30) 

becomes 


The total free energy change on nucleus formation is the sum of 
volume and surface terms 


rdG- 


(30) 


(31) 


ag_ 


iS - 3 IT (T-x^) * 4ir r o 


(32) 


where cr is the liquid-solid interfacial energy for the binary alloy, r 
is the radius of the nucleus, and Vg and are the specific volumes of 
components B and A respectively. Thus the denominator of the first term 
is simply the total system volume. Maximizing AC v/ith respect to r 
leads directly to the critical radius 


Z7 


2oN[X 2 Vg + (l-Xj) v^] 


iG, 


(33) 


m 
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2.1.4 (Continued) 
and hence to the critical free energy 
16110 V[Xj,Vp (1-X„)v„]^ 

aG- (34) 

m 

From equation (34), the nucleation rate in the alloy can now be calculated 
from the usual relation 


I = Iq exp [ 


kT 


] 
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2.2 THEORETICAL TECHNIQUES RELATING TO UNDERCOOLING, EFFECTS 
OF GRAVITY ON MELTS AND SOLIDIFICATION 

2.2.1 The Relation of Undercooling to Intermolecular Potentials 

The undercooling of liquids is of considerable importance in industrial 
solidification processes [16,17], in the study of the liquid-solid 
transition [18-24], and as an interesting phenomena in its ov;n right 
[25-28]. One fundamental question to be answered about the undercooling 
of pure liquids can be stated: "Why does one liquid material undercool 

more (or less) than another?” For simplicity, the following discussion 
of this question will be limited to impurity -free, single-component 
liquids.- 

The most extensive experimental v/ork applicable to the above question is 
the study of maximum undercooling and homogeneous nucleation performed 
twenty-five years ago by Turnbull [23,24], and his data are still the 
most widely quoted in the literature. The Important result of Turnbull's 
experiments, which used the “small drop" technique to ensure purity and 
homogeneous nucleation, was the establishment of a reproducible maximum 
value of undercooling, AT , at which each liquid material nucleates 

niciA 

homogeneously. The question of why these materials nucleate homogeneously 
at a characteristic undercooling has been largely ignored, except for 
discussions of the relation of undercooling and solidification to liquid 
structure [4,25], To state that the amount of undercooling a pure 
material melt will sustain depends on the structure of the solid to 
which it solidifies is. almost equivalent to stating that the undercooling 
depends on the intermolecular potentials in the solid phase of the 
material. It is the purpose of this report to shov; that, within the 
framework of classical homogeneous nucleation theory, the maximum 
undercooling obtainable before a given material solidifies can be 
related (at least formally) <:o the intermolecular potentials in both : 
the liquid and solid phases. ^ 
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RELATION OF UNDEPXOOLIKG TO MACROSCOPIC PROPERTIES 


The classical homogeneous nucleation rate, I 


/ nuclei \ 
^cm^sec^ 




2 


is given by [20] 


( 35 ) 


where, for spherical nuclei. 


A2 = 


16 n cf^ 

c 



(36) 


v/here T is absolute temperature, a is the liquid-solid Interfacial energy, 
S^ is the eni 
Turnbull [4] 


S is the entropy of fusion and k is Boltzmann's constant. According to 


thus 


I = ,n33 = .76 nycM, 
■''o cm^sec i 


I = exp[76 - 


(37) 


Obviously, the largest undercooling a liquid can sustain is the largest 
value of AT which leaves I in equation (37) essentially zero. If one 
takes this value of I to be less than or equal to some small, constant 
value, 5, then from equation (37) 

76 - . 


OfilGINAL PAGK IS 
OJ’ POOR quality 
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2,2.1 (Continued) 


or _ A 


( 38 ) 


Therefore, the larger the value of A, or by equation (36) the larger the 

3/2 

value of , for a given material , the larger undercooling, the 

Op rUaX 

material will sustain, A rigorous test of how well this classical pre- 
diction agrees with experiment must av/ait further application of a direct 
technique for measuring interfacial energies, such as that developed by 
Glicksman and Void [29], This is because most published values of o 
were calculated from nucleatlon data via equation (37) [19], since the 
Glicksman and Void technique has not been applied to many liquids. It 
is no surprise, then, that values of A calculated for various liquids 
using Turnbull's values for a [23] should correlate with measured aT^^ 
values as \vell as they do in Figure 4. Data for calculating A values 
are given in Table I. Note that I and AT,^„ are the variables measured 

by Turnbull, and that S (=r^) is the other “known" parameter used to 

‘m 

calculate a from the nucleatlon experiments. 

RELATION OF TO IMTERMQLECULAR POTENTIALS 


Since entropy is a state function, the entropy of fusion is calculated 
from the partition function by [30] 



“ k In + kT [g^ {In 


(39) 


v/here is the liquid partition function and is the solid partition 
function. In the canonical ensemble, the partition function is given by 
the rigorous relation [30] 



(V,T} = / 




(P.F) 


kT 


(40) 
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TABLE I DATA FOR CALCULATING VALUES OF A (from Turnbull [23]) 


MATERIAL 

0 

erg/cm^ 

Q 

10^ 

cm* 

*^M 

“X 

lf,B St'S 

MB! 

t — 0 

A . 

cm^°K 

Mercury 

28.1 

1.583 

234.6 

6.747 

77 

1:57.8 

612.7 

Gal 1 1 um 

56.0 

4.752 

302,8 

15.693 

76 

226.8 

617,8 

Tin 

59.0 

4.415 

504.9 

8,528 

118 

386.9 

941 .3 

Bismuth 

54.4 

5.107 

544.0 

9,388 

90 

454.0 

698.8 

Lead 

33.3 

2.789 

600.0 

4.340 

80 

520,0 

676.5 

Antimony 

101.0 

: 10.778 

903.0 

11.936 

135 

768.0 

; 1069.0 

A1 umi num 

93.0 

10.669 

931.7 

11.451 

130 

801.7 

963.7 

Germani um 

181.0 

23.38 

1231.7 

18.982 

227 

1004.7 

1410,0. 

Si 1 ver 

126.0 

10.97 

1233.7 

8.892 

227 

1006.7 

1746.6 

Gold 

132.0 

12.429 

1336.0 

9.303 

230 

1106.0 

1707.0 

Copper 

177.0 

18.327 

1356.0 

13;515 

236 

1120.0 

1813.9 

Manganese 

206.0 

19.816 

1493.0 

13.273 

308 

1185.0 

2254.5: 

Nickel 

255.0 

26.689 

. 1725.0 

15.472 

318 

1407.0 

2444.5 

Cobalt 

234.0 

22.731 

1763.0 

12.893 

330 

1433.0 

2555,0 

Iron 

204.0 

21 .627 

1803.0 

11.995 

295 

1508.0 

2179.3 

Palladium 

209.0 

18.804 

1828.0 

10.287 

332 

1496.0 

2645.7 

Platinum 

240.0 

21.61 

2042.0 

10.583 

370 

1672.0 

2993.4 

' 
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2.2.1 (Continued) 

where N is the number of molecules in the solidifying sys tenia h is 
Planck's constant, p represents molecule momenta, r represents molecular 
position and Hamiltonian for the system of molecules 

in the liquid or solid state, respectively. 


In the liquid state, the Hamiltonian is usually written as [31] 


1>J 


(41) 


where ^(r*.) is the potential energy between molecules i and j in the 

i J 

liquid. The intergration over momenta may be performed in a straight- 
forward manner to yield 

Z, = ^ /„ exp [- 1.] d"'’'r 


N ! 


(42) 


where A is the "thermal wavelength" 



The evaluation of the "configuration integral" in equation (42) is the 
central problem of liquid state physics and only approximate solutions 
based on a variety of simplifying assumptions have been obtained to 
date [32], Hov/ever, equations (39) and (42) do show explicitly the 
dependence of on liquid intermolecular potentials. 

A choice of models is also available for evaluation of the partition 
function in the solid state. Perhaps the Einstein model of a solid 
is the simplest model to have been reasonably successful in predicting 
solid properties [33]. In this model , the. partition function is 

^ , 'FT-i 

Zs = TF (l-e ) (43) 
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where the Wi are the phonon frequencies derived from the Hamiltonian for 
the solid lattice of N molecules v/hich is appro.d mated by a sum of terms, 
each representing a harmonic oscillator. The frequences in the harmonic 
oscillator approximation are related to the intermolecular potentials 
by 

=7"'“! >'fj • . w 

This completes the formal relationship between entropy of fusion and 
intermolecular potentials in the solid and liquid states. 

RELATIONSHIP OF a TO INTERMOLECULAR POTENTIALS 

Although the liquid-solid interfacial energy, a, is of major importance 
to most nucleation and solidification models*, it has remained one of 
the most difficult parameters to measure directly by experiment or to 
calculate theoretically. The reasons for the theoretical difficulties 
arise primarily from the inability to describe liquid structure, 
particularly in the interface region. Some recent vfork by Spaepen [34-] 
based on the random packing of hard spheres model developed by Bernal 
and coworkers [35-37], derives an expression for <s in terms of the 

*In common practice, little distinction is made between the interfacial 
energy, a, and the interfacial tension, y ^ - The primary distinction 
between the tv^o is that a is a scalar energy per unit area while y is a 
vector force per unit length. There is also the thermodynamic relationship 

where y| is the magnitude of y. , which further differentiates between 
a and y. 
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configurational entropies of the interface and the bulk material. This 
treatment is strictly structural > and the interface whose properties 
are studied is itself built according to the construction rules which 
apply in the random packed phase. This approach therefore sho'ws no 
explicit relation between o and the intermolecular potentials, but 
certainly these potentials must play some role in establishing the 
random packed structure from v-zhich the construction rules are derived. 

The most obvious formal relationship of a to intermolecular potentials 
is found in the thermodynamic expression for a as the change in free 
energy F, with interface area 




P.T.M 


F is related to the partition function for the system by 


F = -kT In Z 


and hence to the potentials ^ as discussed above. The difficulty occurs 
in deriving F{A) for a system containing a liquid-solid interface. Thus 
alternate approaches to the derivation of the liquid-solid interfacial 
energy have been pursued historically. 

The theory which shows the most direct relationship between liquid-solid 
interfacial energy and intevTiiDlecular potentials is that due to Skapski 
[38]. In this theory the difference between intermolecular potentials 
in the solid and liquid phases appears explicitly ip the equation for a. 
In actual computations, the difference in intermolecular potentials Is 
replaced by a function of the heat of fusion. Skapski 's theory, although 
simple, yields values of 0 v/hich are in order-of-magnitude agreement with 
those obtained by Turnbull from nucTeation experiments (see Table I). 
Hov/ever, both this theory and that developed by Turnbull [39], which is 
similar to Skapski 's theory in the dependence of 0 on intermolecular 
potentials, are based on a general randomness of the liquid (and perfect 
crystallinity of the solid), and use no specific structural information. 
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2.2.1 (Continued) 

A n'gov'ous analysis of liquid-solid interfacial energy given by Herring 
[40] was used by Glicksman and Void [29] as a basis for their method of 
experimentally determining cr. This method is of particular importance 
since it is at present the only source of experimental values of a 
which are not obtained from nucleation experiments. The theory, as 
applied by Glicksman and Void, considers a symmetrical low-angle tilt 
boundary intersecting the liquid-solid interface under isothermal, stress- 
free conditions. In this case, the tilt boundary must terminate at a 
symmetrical groove on the liquid-solid interface. If f is the 
dihedral angle subtended by the liquid phase at the root of the groove, 
Glicksman and Void show that is related to the tilt angle, e, and the 
liquid-solid interfacial energy by 

1 I E 

j cos = A" - '^lyi 0 (45) 


where A" is a constant. When the angles ijt and e are measured experimentally, 
the parameter — cos may be plotted as a linear function of In 0 with a 
slope 


Thus 



(46) 


where m is determined empirically. is a factor of the Read-Shockley 
formulation for the internal energy per unit area of a symmetric low- 
angle tilt boundary, and is given by 


c — G a 
0 4it(T-v) 


(cos e + sin e) 


(47) 


v/here e is the orientation angle for the boundary in a selected coordinate 
system, a is the lattice constant, v is Poisson's ratio and G is the 


2-24 



D256-10202 


2.2.1 (Continued) 

modulus of rigidity. It should be noted here that in their theory of 
crystal-melt interfacial energy, Kotze and Kuhlmann-Wilsdorf [41] derive 
the relation 

o = 0.85 [-J-+ 2] 

m loir 1-v (48 

where b is the Burger's vector and is a constant angle, approximately 
25°, The derivation of this equation considers the liquid-solid interface 
to be half of a general dislocation grain boundary, a central idea in 
the Kuhlmann-Vlilsdorf theory of melting [42]. 

The importance of equation 48 and of equations 46 and 47 together is the 
direct dependance of a on the modulus of rigidity in the solid, 6. In 
most real materials, anisotropy cannot be neglected. Therefore, G 
will depend on the elastic constants of the material. In the Kuhlmann- 
Wilsdorf model the relationship is 



3 C 


C]! - 


^ ^44 " ^12 


(49) 


while the relationship developed by Read and Shockley [43] whose theory 
is used by Glicksman and Void is 


G 

1-v 



+ C 


12 


^ r 44 /^ 11 : “ ^12 
C]! 


)] 


T/2 


(50) 
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Tlie elastic constants, are the traditional form, obtained by direct 

contraction, of the tensor elastic constants, C _ . , It is the C „ , 

aP,YA a0,YA 

which are directly related to the intennolecular potentials, being 
complicated linear combinations of second derivations of the configurational 
energy function 

® (51) 

in which <{t(|K.-x.|) is the to-body potential betv/een molecules at position 

• J 

X. and at position x-. According to Huntington [44], the dependence of 

1 J 

C „ , on ® takes the form 
aS,Y^ 

‘^aY,3A ^aA,0Y " ^ («Y»PX) + (aX,3Y) (52) 

where 

[ap,YA] 
and 

(aY>l3A) 


Here v represents the volume per unit cell and the x are the components 
of the relative position vector between the lattice positions associated 
v/ith the indices given in the curly brackets. is the square 

matrix which is inverse to the matrix 

V r-f f 

p 3X 3X k‘'^ 

f. u V 
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Thus the Interfacial energy is seen to ultimately depend on the second 
derivative of the interniolecular potentials in tMr solid phase. Knov/ledge 
of what influence the liquid phase potentials may have on the interface 
structure and interfacial energy was precluded by the formulation of the 
Herring theory in terms of an equilibrium betv;een the solid and liquid 
which leads to zero difference in free energy across the interface. In 
the Kuhlmann’-Wilsdorf theory, the liquid is viewed as a solid with a high 
concentration of dislocation cores, with the interface assumed to be a 
modified dislocation boundary, hence the obvious dependence on solid 
properties. What is needed for an understanding of the fundamental 
origins of liquid-solid Interfacial energy is a theory which combines 
the liquid structural considerations of Spaepen with the solid potential 
energy relationships demonstrated above. 

SUMMARY 

It has been shown, within the framework of classical physics (with the 
possible exception of the Einstein model), that the amount of undercooling 
the liquid phase of a pure material can sustain is a function of the 
intermolecLilar potentials of the material. That this conclusion may also 
be valid for heterogeneous nucleation on impurities or substrates is 
indicated by the postulated dependence of undercooling on the disreglstry 
[12,16,45]. It might be noted that Burton [46] and others [8] have 
applied statistical mechanical considerations similar to those above to 
the general problem of nucleation from the vapor, but have not addressed 
the liquid-solid transition. The present work illustrates the difficulty 
In carrying out rigorous calculations of the characteristic maximum 
undercooling to be expected in real liquid systems based only on knov/ledge 
of Intermolecular potentials, without making simplifying assumptions similar 
to those of Skapski in his nsarest-neighbors model of the liquid-solid 
interface. 
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Although mors realistic nucleation models than classical homogeneous 
nucleation theory could be used in determining characteristic maximum 
undercooling, predictions must eventually depend on thermodynamic values 
like cj and S . Until these quantities can be rigorously evaluated from 
intermolecular potentials for various materials, the formal relationships 
presented here must suffice as the explanation for v/hy different pure 
liquids undercool different amounts. 

2.2.2 Undercooling and Growth Rate 


That the growth rate of dendrites and other solid interfaces into under- 
cooled melts is proportional to the square of the undercooling is well 
established experimentally [47]. Theoretically, there are several ways 
to show the relationship. An early method, due to Chalmers [48], addresses 
the problem of lamellar growth of solid into the melt. By deriving re- 
lationship between geometric, kinetic and energetic properties of the 
lamellae, Chalmers was able to show that 



(55) 


where is growth rate, K 
{■|r = one lamellar width) > 

>-‘sr 


^ is a constant and X is the lamellar wavelength 
which is related to the undercooling by [4s] 

(56) 


Thus 



wi th 




K. 


8 D S 


K, " crnrrc 


a 


E 


(57) 


(58) 
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where D = diffusion coefficient 

S = entropy of fusion 
- liquidus slope 
k = distribution coefficient 

a 

= equilibrium concentration of the melt 
" interfacial energy at the a-3 lamella phase boundary. 

The more general solidification rate theory of Turnbull, Hillig [22] 
and Cohen [49] derives the expression 

(59} 


(60) 


( 61 ) 


and other terms are as defined above. These simple relationships between 
growth rate and undercooling are complicated by the lack of good data 
for the parameters which make up K^/K^ and A^/A-j , but are nevertheless 
most useful in the analysis and guidance of solidification experiments, 
as v.'ill be demonstrated in Section 2.3. 
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Uc = ^ [l-e“^2AT^ 

Expansion of the exponential and multiplication yields 


plus higher order terms, with 

A„ V D 
2 „ c 

A-j 4 'st kT a 

where v = specific volume 

k “ Boltzmann's constant 
o = liquid-solid interfacial energy 
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2.2.3 Interface Stability and Dendrite Formation 

The way undercooling influences solid structure during growth is 
illustrated by Figure 5. 



FIGURE 5 


THERMAL GRADIENTS AT INTERFACE 
(from Reference 16) 


When the melt is undercooled, the thermal gradient in front of the interface 




X 


X, 


m 


(62) 


is negative since T < T^^. A negative thermal gradient implies a negative 
gradient of free energy in this region. If an instability occurs in the 
interface, i,e., the formation of a protuberance, the protuberance vjould 
grow into the region of negative vT whereas it would merely melt if ^T 
were positive. A protuberance which grows becomes a dendrite, and as 
long as undercooling is maintained in the melt,' dendritic growth is 
sustained. In the case of eutectic alloys, an interface instability also 
results in lamellar growth. These dendritic or lamellar structures are 
highly anisotropic, a result which may be desirable or detrimental, 
depending on the desired application of the resultant solid. 
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A fornial theory of liquid-solid interface stability has been developed 
by Sekerka [50,51] and extended by Chen [52] to include semi-static 
external field effects. The geometry of an unstable interface is 
illustrated by Figure 6 - the dotted line representing a stable, planar 
interface moving into the liquid with assumed constant velocity V, and 
the curved line representing a possible shape of an unstable interface 
in the process of forming a cellular structure in the solid v/ith cells 
of compositions a and b forming behind the maxima and minima, respectively, 
of the unstable interface. 

It is assumed that the interface is initially planar, and that at some 
time a small perturbation <f!(x,t) is introduced Into the process such that 

z ~ i|)Cx,t} = 6(t) sin w X (63) 

where u - 2tt/X with X being the wavelength of the perturbation, or the 
distance over which the concentration varies on a microscopic scale. 
Obviously, if the perturbation amplitude, 6', increases with time (6>0), 
the perturbation will grow, the interface will be unstable and the 
cellular structure, or micY'‘osegregation, will result. But if 6 does 
not increase with time (5£0), the interface will remain stable and no 
microsegregation will occur. So the problem becomes the determination 
of 6 as a function of the solidification parameters, in particular, the 
solute molecule velocity u and the cell size X (or u) . u is presumed 
to be induced by an external field, F, rather than being due to any in- 
ternal process, such as Brownian motion. 

Calculation of 6 is begun by solving the steady state solidification 
equations 

v'c+ff -0 (64) 


=0 

3Z 


(65) 
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FIGURE 6 INTERFACE GEOMETRY FOR SEKERKA THEORY,' 


T • I l V 8T’ _ n 

in sol id V T' -{* — -rr — 0 

“s 

subject to the boundary conditions 


( 66 ) 




v(x) = V + 5 sin ux = 

1 [K, (|f) - K, (|1)^] 
P ^ 1 


(67) 


(68) 


where C is solute concentration V ci is thermal diffusivity, m is the slope 
of the liquidus curve for the system (assumed negative) » T^^ is the melting 
temperrvtures Q is the solvent latent heat of fusion,: a is the intcrfacial 
free iivorgy, K is the Interface curvature, k is the segregation coefficient 
(0<k = Cg/C|^<l) and V‘ = V*-u, The left portion of equation (68) states 
that the real interface or solidification velocity, v(x) is identically 
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the sum of the macroscopic planar interface velocitys Vj, and the time 
derivative of the small variation fron, planar geometry, ^(x,t). The 
right portion of (68) states that the solidification velocity at each 
point X on the interface when calculated from heat flow considerations 
must agree with the velocity calculated from diffusion considerations. 

The solution for 5 obtained from equations (64) -{68} is an extremely 
complicated function of the parameters listed above as v/ell as the 
temperature gradients (which are assumed positive) and thermal con- 
ductivities in the liquid and solid, the concentration gradient and 
diffusion coefficient in the liquid and the parameter (V-u/p) where 
p = 1-k. The sign of 6, and thus the occurrence of microsegregation, 
will depend on the sign of the expression vjhich is obtained for 6/6 

5 / 6 = ( 69 ) 

where the denominator, D(tj) and the numerator, N(o)) both depend on 
(V-u/p) as well as w and other "fixed” solidification parameters. There- 
fore, the solute velocity, which can depend on external fields, enters 
the determining expression for microsegregation through the term (V-u/p), 
thus defining the three cases of interest: 

V - - = 0 

P 

V - -JJ- < 0 

p 

V-f>0 . 

The analysis which has been performed [52] finds that the interface is 
stable (no microsegregation) for V-“_< 0. When V-~ > 0, the value of 
6 depends on the relationship of V- to the parameter 

k T rV y 
R - m . 

(-m) DpC^ 
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where r = c/Q and where ts liquid concentration in the flat interface 
case. If V-^ 5 B, the interface remains stable. When V-— > Bj the 
stability criteria are very complicated [52] > but correspond to the 
stability analysis given by Sekerka [51] if the parameters k and V are 
redefined as k' = k V/V and V’ = V-u. Thus for negative m, positive 
temperature gradients and for 0<k<l, microsegregation of solute in the 
solid cannot occur for V-u/p <. B, and for microsegregation to occur 
when V-u/p > B, the system must satisfy the instability criteria described 
by Sekerka [51]. 

There remains the question of segregation cell size and its relation to 
the occurrence of microsegregation. If the entire equation for 6/6 is 
plotted versus for a given value of V-u/p, curves such as those in 
Figure 7 result. Curve 1, in which microsegregation is seen to occur, 



FIGURE 7 SCHEMATIC CURVES OF 6/6 

indicates that the cell size will Tie between “ 1 ^* 

When the tedious calculation is carried out, it is found that the range 

-3 

from ^2 ^-1 order of 5 X 10 cm, v/hich is in reasonable 

agreement with experiment [50]. Beyond this order^of -magnitude 
estimation, further accuracy is impossible because Sekerka 's theory 
is based on (a) steady-state approximations, and (b) assumed small 
amplitude (linear) perturbations. A transient theory has been 
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developed [53] to eliminate inaccuracies due to (a), but because of the 
assumptions in {b)a the transient theory does not yield any significant 
improvement in accuracy. Thus it seems that the indication of the 
occurrence of microsegregation by the sign of 6, and the approach to cell 
size determination exemplified by Figure 7 constitute the best results 
available from a theoretical model v/hich describes interface stability 
and predicts solute distribution characteristics in solids from liquid 
state parameters. 

2.2.4 Gravitational Settling of Nuclei 

When a multiphase alloy solidifies in a gravitational field, the 
settling of nuclei of the heavier phase(s) toward the bottom of the melt 
results in macroscopic segregation of this solute phase in the solid 
material. The distribution of nuclei of the heavier phase, C, in the 
melt should be described by Pick's law 

D ^ = - CU (70) 

where D is a "nucleus diffusion coefficient", z is the direction of the 
gravity vector and also the axial direction in an assumed cyltndrically 
shaped melt and U is nucleus settling velocity. The equation of motion 
for the nucleus, including friction and buoyant forces is thus 

H ^ = Mg - SwirU - |r r^(z)g " (71) 

where li is nucleus mass, r its radius, n the melt viscosity and p(z) the 

melt density (which is assumed to change more si ov/ly with time than U), 

2 

g is the acceleration due to gravity, 930 cni/sec . When equation 71 is 
solved for U, the result may be written as a function of p(z) and g 

U(z,t) = [1 ygD(l-e“^/^'^) (72) 
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where pj^ = snd y - . For most reasonable nucleus motion, 

t » yD, so the exponential vanishes, leaving (for slowly varying 
melt density) 

U(z) = Y g D [1 - 4^] . (73) 

Pm 

Substitution of this result into equation 70 and performance of the 
integration wi ‘;h respect to z leads to the integral equation 


InC = - Y g 2 + y g f (z) 


where 


p[C(y)] = 


P [C(y)] dy and 

0 

nil C(y) - iHg C{y) t 
v-j C(y) - Vg C(y) + V 2 C.J- 


(74) 


(75) 


with subscripts 1 and 2 referring to solute and solvent atoms, respectively, 
and Cj being the total (homogeneous) melt concentration or nifmber density. 

By combining constants, this can be written 


A. C + tn' 

p(C) = 

and by noting that equation 74 can be written 

C = exp [yg (f-z)] , 

clC cfi* 

one can calculate ^ and solve for - 

df = 1 + J_ dC_ 
di Y 9 C dz 


(76) 


(77) 
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ci 'p 

But from the definition of f(z), is also 

^ - 1 ^3 ^ 

dz " pj,j C + V' 


( 78 ) 


Equating the right-hand sides of (77) and (78) leads to the new differential 
equation ■ 


(A. C + v') dC 

— 2" — = YSZ 

BgC 


with 



and 




(79) 


This can be integrated directly, and applying the boundary condition 
that C = (the solute concentration equals the homogeneous concentration) 
at 2 = 0 (the center plane of the melt), the solution obtained is 


L^-^l 

B C + ^1 h 

( . } 4 -) 


(™) 


v' 


eY9^ 


(80) 


This is a complicated function of both concentration and melt physical 
parameters. Therefore its behavior is best illustrated by application to 
specific problems where all the constants on the left-hand side can be 
evaluated, as will be done in Section 2.3.2. It should be noted, hov/ever, 
that since the parameter y depends on the mass and thus the size of the 
nuclei whiclvare settling, it therefore depends on (^.T) , and adds further 
complexity to equation 80. 
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2.2.5 D'irect (Non-Convection) Effects of External Fields on Melts 

It has been shown [18] that when a liquid is subjected to changes in 
external fields acting on it, the therinodynami c variables most directly 
affected are volume and temperature. In parti cularj if the gravitational 
field acting on a melt changes by an amount fig, a change in hydrostatic 
pressure in the melt will result, inducing a corresponding change in 
specific volume, v, given by 

= 3 p z 6g (81) 

where a is the isothermal compressibility, p is the density of the melt> 
and z is the distance below the liquid free surface. The pressure 
change also induces a temperature change given by 

^ f P ^ S3 (82) 

where a is the thermal expansion coefficient. Other variables affected 
by a field change through their dependence on volume or temperature 
include: diffusion coefficient, Gibbs free energy, and interfacial 

energy, if a is dependent on T in the first place. More complicated 
parameters, such as nucleation rate, nucleus radius, solidification rate, 
total solidification time and grain size depend on combinations of the 
above variables. The heterogeneous nucleation rate is one of the most 
important parameters affected by field changes, and the calculation of the 
altered nucleation rate, I', is greatly simplified if the interfacial 
energies are not influenced by field changes - that is, if a does not 
depend strongly on temperature. The problem can be illustrated by 
recalling that I is an exponential function of f(e) where 

f(o) = j (1-cos 9)^ (2-fcos e) 
and 

cos 
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where the subscripts s, I and n refer to substrate, liquid and 
nucleus, respectively. If a., is a function of T, and therefore, 

^ tJ 

under changed field conditions, of 6T, the difficulty of calculating 
I'(<ST) over calculating I' for constant f(e) or fixed a... Is obvious 

* u 

At present there is no experimental evidence and no theoretical 
justification, other than Spaepen's work [34], for considering a or 
o.. to be more than very weakly dependent on temperature. 

I J 

In general, the effect which external field changes have on physical 
properties of melts, such as those listed above, may be found by 
replacing T and v in the equations for the properties of interest 
with T' and u', whev'o the primes V'efer to the values of the variables 
in the changed field: 

r = T t 5T 


and 

u ' = y + 6u . 

The field effect on undercooling is thus 
AT' = Tj^ - T' == T^,^ - T - 6T - aT-6T . 


Hence the undercooling is decreased by a positive 6T and enhanced 
by a negative 6T. It should be noted that free energy calculations 
require more knowledge of the field conditions than is contained in 
equations 83-85, due to the direct dependence on pressure. Also, 
equations 81 and 82, as they ai'e written, hold only for gravity field 
changes. For iiiagnetic fields the corresponding equations are [18] 


ov 

V 


1 

2 








2-39 



D256-10202 


2.2.5 (Continued) 
and 

v/here is the permeability of free space, x is magnetic susceptibility, 
is the initial and H the final field strength in the liquid, and Q 
is heat of fusion. 

2.2.6 The Relation of Solid Physical Properties to Structure and 
Electron and Phonon Scattering 

The influence of undercooling, nucleation and general solidification 
history on the physical properties of materials solidifying from the 
melt was investigated as part of the attempt to understand how all 
solid properties might be better predicted and controlled by manipulation 
of the solidification process. It appears that no direct link exists 
betv/een the solidification process and the physical properties of the 
resultant solid. Hov/ever, the two are linked throOgh the structural nature 
of the solid which depends directly on solidification history, primarily 
the undercooling and nucleation of the melt. Many physical properties 
of solids, such as thermal conductivity, electrical resistance, Grilneisen 
parameter, and specific heat depend on the scattering of electrons and 
phonons by non-periodic structures in the solid. In particular, three 
physical properties of solids, the thermal conductivity, k, the 
electrical conductivity, cr, and the thermoelectric pov-/er, Q‘, have been 
found to depend on the phonon mean free path, A, in the solid. Specifically 
the relationships are given by Ziman [54,55] as: 

( 88 ) 

(89) 

(90) 

3 ■ 


K = y V C A 


a “ 


3e^ 

TT k 1 


K, 


Q' y.) 
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vn‘th 1<2 and being the magnitudes of the tensor transport coefficients 




/A 


-> ~y 
V V 

e u V 



(91) 


(92) 

= 1 Tr^(kT)^ ?Q (u) ' (93) 


where v 
c 
e 



= average phonon velocity in the solid 
- heat capacity 
= electronic charge 
= differential of the Fermi surface 
= absolute temperature 
= velocity of the electron or hole carrier 
= I ? I 

= electron mean free path in the solid 


and y is the free energy of an electrons i.e., its chemical potential. 
Ziman also reports [55] that the phonon mean free path in polycrystalline 
solids is on the order of the size of the grains. That is, 



,3 / 

/V 

V g 


(94) 


th 

where v^ is the volume of the i crystal grain. Thus grain size in a 
polycrystalline solid can be expected to have some effect on the thermal 
conductivity of the solid, and through k, on a and Q' as well. This 
effect is usually masked by the considerable temperature dependence of 
K through c. As is well known, structural properties such as strength, 
hardness and plastic deformation are more directly related to solid 
structure and hence to the solidification process. The relationship 
of the solidification environmeni: and process to structural and physical 


properties is illustrated by Figure 8. 
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2,2.6 (Continued) 



FIGURE 8 HYPOTHETICAL RELATION OF SOLIDIFICATIOri PROCESSES 
TO PHYSICAL PROPERTIES 

2.2.7 A Solidification Model to Predict Grain Size and Dendrite 
Morphology 

For a long time, it has been known empirically that grain size in a poly- 
crystalline solid depended on solidification history, primarily on under- 
cooling. More recently, it was shown experimentally [56] that the 
morphology of dendrites also depends on these processes. If the proper 
mathematical relationships betVN'een nucleation and growth of crystalline 
grains from the melt can be found, and if the proper external boundary 
conditions can be applied, then in principle it should be possible to 
predict something of a solid's structure from a knowledge of the melt 
properties and environment. An attempt to develop such a predictive 
model will now be described. 
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2.2.7 (Continued) 

Kattamis and Flemings [56] reported observing a distinct change in 
dendrite morphology as a function of melt undercooling. To describe 
this theoretically, a "morphology parameter" may be defined, which 
is simply the ratio of that portion of a dendrite's volume due to 
grov;th, v^, to the portion which nucleates initially, v^^. Figure 9 
shows a cylindrical dendrite growing from a spherical cap nucleus formed 
on a heterogeneous nucleant substrate, characterized by v/etting angle e. 



FIGURE 9 SIMPLE DEMDRITE - NUCLEUS MODEL 

As was shown in Section 2.1,3, the nucleus volume is given by 

v^=|Trr\{G) . (95) 

The area of the dendrite base is that of the nucleus and is simply 
V r sin G. To a close approximation, the dendrite length will be the 
product of the growth rate, with the growing time, dt^, which is 
the time required for solute atoms to diffuse from the extremities 
of the region in which a dendrite is growing to the tip of the dendrite. 
Hence 

M = T, r^sin^o U 6t , . (96) 

9 c u 
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2,2.7 (Continued) 

Mow r depends on undercooling via equation 19, f(6) and sin^e depend 
on the nucleant, l)^ is a function of undercooling through equations 
60 and 61, and 6t^j is given by 

••d -/r ij>'" <”> 

where ^ is the volume of a grain and 
71 in p a^ 

with in = solute atom mass 

= solute atom density 
a = solute atom radius 
k “ Boltzmann's constant. 


The total volume of material, V, is assumed known, so the remaining 
unknown is N,the total number of nuclei which form in the melt. This 
number is 


M = r V 5t^ 


(98a) 


or the grain size is 


i = 

W I6tj^ 


(98b) 


where I is nucleation rate (equations 1 and 20) and fit is the time 
required for heat to flow' out of that region of the melt where a single 
nucleus is forming and is expressed by 





(99) 


where is the thermal diffusivity of the melt. In this first attempt 
at development of a predictive solidification model , equations 97 and 99, 
the diffusion and heat flow constraints, are derived by dividing the melt 
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2.2.7 (Continued) 

into N imaginary cubes. One nucleus is assumed to form at the center 
of each cube at some time during the interval 5tj^ > st^, which is taken 
as the total solidification time. The maximum distance which either a 
solute molecule or heat must travel is that distance from the center of 
the cube to a corner, ^ . Dimensional analysis of the fundamental 

diffusion and heat flow equations then leads to equations 97 and 99, 

Equations 95 through 99 plus 1, 19, 20, 60 and 61 may be solved simul- 
taneously (a Wang 700A computer program has been developed for this 
purpose) to yield values of all unknown parameters, particularly the 
morphology, v^/v^, and the grain size, ^ . Of course ^is really the 
size of each imaginary cube, but its value gives a reasonable indication 
of average grain size in the solid and is an order of magnitude 

estimate of dendrite spacing. There are three value ranges of v /v 

g n 

which are of interest, each corresponding to a morphology observed 
by Kattamis and Flemings, and shown in Table II. 


TABLE II DENDRITE MORPHOLOGY CORRELATION 


Value of Vg/Vj^ 

Dendrite Morphology 

» 1 

normal 


cylindrical 

« 1 

spherical 


This crude first cu’ at a mathematical model to calculate grain size 
and morphology in solids from solidification history contains several 
simplifying assumptions or deficiencies which should be remedied by 
more rigorous analyses. The most obvious of these are: 

1, No accoLint is taken of dendrite branching (other than as a heterogeneous 
nucleation site for a new dendrite) or of the limiting of one dendrite's 
growth by other dendrites (except to limit each dendrite to its cube). 
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2.2.7 (Continued) 

2. The heat flov; constraint treats only the local (cube) thermal 
environment, and treats each cube independently; total melt heat 
transfer analysis is required. 

3. A better definition of total solidification time is required. 

4. Model accuracy is no better than the “constant" input parameters 
used for numerical calculation. At present, very few of the 
properties of melts required by this model are known v/ith any 
accuracy. 

Although these obstacles are formidible, calculations can be made v;ith 
the present model which are accurate enough to at least allov/ comparison 
with experiment, as is done in Section 2.3.3. The reason for pursuing 
development of this model is clear: it forms the core of a much larger 

model v/hose purpose is nothing less than complete analytical materials 
design for user requirements, using gravity, or microgravity, as a 
primary variable. Such a model will be discussed in Section 2.4.1. 
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2.3 APPLICATION OF THEORY TO EXPERIMENTAL PROBLEMS 

2.3.1 Undercooling and Growth of NH^CI-HgO in Microgravity 

Two solidification experiments vn'th the metal -model material NH^Cl-H^O 
(ammonium chloride solution) have been flov/n by Johnston and Griner of 
MSEC in the Space Processing Applications Rocket (SPAR) experiment program 
The December, 1975 flight experiment sample began to solidify at a knovm 
grov/th rate, and it v/?s useful to calculate the undercooling corresponding 
to that growth rate. The May, 1976 flight experiment sample undercooled 
a known amount, but failed to solidify. Appropriate portions of the 
theory from Sections 2,1 and 2.2 were applied to these experiments to 
provide further insight into vjhat occurred experimentally. 

In order to perform calculations of nucleation and growth rates, it is 
necessary to have values of melt parameters such as entropy of fusion, 

S , liquid-solid interfacial energy, a, and diffusion coefficient, D. 
Unfortunately, none of these values is known for NH^iCl-H20. However, 
being a metal -model material, it vnll have values of S^, a and D which 
are approximately in the range common to most metals. In particular, 
most metals have a value of which is, within about 50?i, equal to 
Boltzmann's constant, 8.3167 X 10^ I 

empirically that a reasonable approximation for cr can be derived from 
the value of S by the relation 

° = „l/3 ,2/3 

where T^^ is the equilibrium melting temparatiire, N is Avogadro's number 
and.v is the specific volume. For the ammonium chloride solution used 
(28“/ by w'eight NH^Cl), the known parameters are 

1. = 295.8 

m- ■ ■ ■ ■ ‘ 

3 

v 6.68 cm /mole . 

2 

Using this information, the Interfacial energy v/as found to be 37 erg/cm . 

-5 2 

Since the diffusion coefficient for most metals is near 4 X 10 cm /sec, 
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2.3.1 (Continued) 

this leads to a value for A«/At in equations 60 and 61 of 0.0036 cm/sec'^K^. 

d, I 

If one writes the heterogeneous nucleation rate equation (see equations 
1 and 20) as 

2 

I = exp [75.985 L liP . L ] (TOOa) 

(T -at) 

' m 

2 

the value of A obtained with the above values of the NH^CI-H^O parameters 
is 39.665 X 10^ 

December 1975 SPAR Experiment 

Since the growth rate was observed to be 0.018 cm/sec (1.08 cm/niin) in 
this experl man tj the problem was to determine the probable amount of 
undercooling causing this rate. This was done by using equation 60 

U = (.0036 aT^ 

see 

to calculate the curve shown in Figure 10. As seen from this curve > 

the undercooling value corresponding to 1.03 cm/min growth rate is 

2.25°K. Metal values of D and S„ do not vary much from those assumed 

2 c 

above, so the product D will probably not affect the calculated 

^ 2 
value of undercooling. Hov/ever, a can vary from about 25 erg/cm to 

30D erg/cm , which leads to the following spread In calculated under- 
cooling values for NH^C 1 “H 20 : 

1°K < at < rK - 

Thus the undercooling, though quite small, was sufficient to initiate 
solidification In the sample in the time available on the rocket flight. 

May 1975 SPAR Experiment 

In this experiment, good temperature versus time flight data curves 
were obtained, but the sample failed to solidify. For the purpose of 
calculation, time data read from these curves is presented (for each 
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AT - °K 

FIGURE 10 GROWTH RATE VS. UflDERCOOLING FOR AMMQMIUM CHLORIDE SOLUTTOi\’ 
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thermistor) in Table III as the "time during which the sample region 

measured by a given thermistor was below the value of undercooling 

shown in column 3". Since the sample ms approximately 5 cm long by 
2 

1 cm ^ a minimum allowable growth rate to solidify in the available 
time can be estimated. The growth rate calculated from the undercooling 
is shown in column 5 , and the time required to solidify the sample at 
that growth rate is in column 6. This data indicates that there was 


TABLE III may ’76 SPAR NH^Cl EXPERIMENT DATA 


Thermistor 

Time 

Sec. 

aT 

Min. U 
cm/sec° 

cm/sec 

Time to Solidify 
sec 

1 

250 

6.8 

.02 

.167 

30 

2 

sample 

heated 

— 

— 


— 

3 

200 

5.8 

.025 

.122 

41 

4 & 5 

240 

17.8 

.021 

1.146 

4.4 

6 

sample 

heated 

— 

--- 

— 

— 

7 

195 

7.8 

.0256 

. "cc ■ 

23 


sufficient time for solidification rf nucleation had occurred. The fact 
that nucleation did not occur at the undercooling measured by the thermistors > 
while it did occur at lower undercooling in the December experiment indicates 
that certain "low angle" heterogeneous nucleating agents v;ere destroyed 
in the melt prior to launch, i.e., the sample was "denucleated" , This 
statement can be quantified if equation 100a is solved for f(e), 

f(o) = [75.9S5 - In I] Al . 

39-665 X 10° 

values of nucleation rate are assumed, and a table of values of 6 versus 
f{fe) is available. Turnbull [7] defines a "measurable" nucleation rate 
as ' 

1 nucUgs ^ j ^ ,t)6 . nu clei ^ 
cia sec cm sec 
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2.3.1 (Continued) 

To assure not only measurable nucleation rates but solidification of the 

9 0 3 

sample within the time available, values for I of 10 to 10 nuclei /cm sec 
v^fere substituted in equation 100b to produce the curves in Figure 11. 

For such a large variation in the band separating the two regions 
indicated is quite small. Since nucleation did not occur in this 
experiment, and undercoolings of up to 17.8°K were obtained in some 
regions of the sample, it is clear that nucleants characterized by 
wetting angles of up to 60° were dissolved into the melt, i.e,, de- 
nucleated. As used here, the wetting angle, e, is a convenient parameter 
for characterizing heterogeneous nucleants, and is not necessarily a 
strict physical quantity. If this were so, every nucleus that formed 
would have to be a spherical cap on a planar impurity particle substrate. 
While this probably occurs some of the time, assuming that all nuclei 
form this way is unjustified. 

That all nucleants v^ith wetting angles lower than 60° were denucleated 
is indicated by the fact that if a nucleus had formed in the melt region 
where aT = 17.8°K, its growth rate, although slower in other parts of 
the melt as shov/n in Table III, would still have been sufficient to 
solidify the sample. That nucleants of angles up to 50° were denucleated 
only in that one highly undercooled region but still existed elsewhere 
in the melt is most unlikely. If nucleants characterized by wetting 
angles of around 62° or more were still active in the melt, a sustained 
undercooling of more than 21 °K would have been required, according 
to Figure 11. The melt regions around thermistors 4 and 5 attained this 
undercooling for about tv/o minutes, indicating that even these nucleants 
had been denucleated. 

2.3.2 Undercooling and Nucleation of Lead-Tin Eutectic Haterial in 
a Centrifuge 

Effects of undercooling and gravity on solidifying materials were studied 
in a series of experiments on Pb-Sn eutectic material by Johnston and 
Griner at MSFC [57]. This material was selected for the drastic difference 


2-51 



WETTING ANGLE 


D256-10202 


2.3.2 (Continued) 


in density of its components. In the experiments, initially solid samples 

of eutectic composition (X == 26 atomic % Pb) were melted in a 1.5 cm long 

by 1.2 cm diameter cylindrical crucible, and held in the molten state v;hile 

inside a centrifuge providing a fixed centripetal acceleration equal to 

2 

some integral multiple of g, 980 cm/sec . The samples were then cooled, 

solidified, removed from the crucible and analyzed metal lographi cal ly to 

determine the volume fraction of Pb-rich dendrites at the bottom of the 

sample. Table IV presents the results of these experiments, along with 

the conversion of the volume fraction, V , to total lead atomic fraction 

* D D 

through the Pb atomic fraction in the dendrites, X . X ms obtained from 
an approximation of the Pb-Sn phase diagram in the temperature-concentration 
region of interest, shown in Figure 12. The value of component concentration 
is related to atomic fraction by 


X = C/C^ 


( 101 ) 


v;here Cy is total melt concentration. Any concentration value is computed 
from 


v.'here Nq is Avogadro's number, 
is its atomic weight. For the 
v/ere calculated 


( 102 ) 

is density of the i^^ component and 
Pb-Sn eutectic malt, the following values 


Cj = 0.35802 X 10^^ cm“^ 

= 0.1038 X 10^^ cm-^ . 

0 


After obtaining other melt parameters, the constants in equation SO can 
be evaUfated as follov/s 


v' = 

m‘ - ^'^Sn^T ^ gm/cm^ 
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TABLE IV HSFC, CENTRIFUGE EXPERIMENT DATA 


SAMPLE 

g 

LEVEL 

AT, 

. 

% 

Pb CONCENTRATION 
IN DENDRITES, XD 

TOTAL Pb 
CONCENTRATION 
ATOMIC % 

2A1 

1 

0 

27 

.710 

38 

3/^1 0 


7.5 

35 

1 

.733 

43 

5A1 

I 

8.6 

37 

-737 

44 

lAl 

1 

n.34 

40.5 

.745 

: 46 

IA3 


0 

46-5 

.710 

47 

2A3 

3 

6-05 

49.5 

.728 

49 

— A 

3 

10.8 

— 

.743 

— 

3A3 

3 

11.34 

56.5 

.745 

53 

5A3 

3 

14.7 

59 

.755 

55 

4A5 

5 

1 

0 

50 

.710 

48.5 

3A5 

5 

1.89 

— 

.715 


2A5 p- 

5 

6.05 

55.5 i 

.728 

52 

1A5 L-J 

5 

7.0 

— 

.762 

{ 

1 

5A5 

5 

8.69 

60 

.737 

1 55 

7A5 

5 

16.44 

45.5 

.760 

! 

I 
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2,3.2 (Continued) 

\ " 0*335 X 10”^^ cm^ 

-23 

^3 ” ^Pb " ^Sn “ ^^*59 X 10 gm 

= ni7p[^ ~ v' = -0.3442 

“ V^N ” \ ^ 

vVBg ^ “2.8 
A^/B^ =0.35 


Equation 80 now becomes 

r(.9577 X IQ-^^) C - .3442^^-^^^ ^1038 X ] 0^3 2.803 _ 

^ - 0.24479 ^ ^ C ^ ® 


(103) 


Since the parameter y depends on the mass and thus the size of the nuclei 
which are being segregated, it therefore depends on the cube of the 
undercool ing» AT. Thus to calculate the lead distribution in a sample 
of size -.75 cm < z < + .75 cm (as shown in Figure 13), particular 



Figure 13 SAHPLE COORDIfJATE 
SYSTEM 


values of both AT and g must be 
specified. If a value of 10®K is 
chosen for AT, curves for lead 
concentration distribution for different 
gravity levels, as shov/n in Figure 14, 
are generated by equation 103. Since 
concentration was measured at only 
one position (z = -0,75) for the MS PC 
samples, no comparison betvjeen theory 
and the MSFC experiment can be made. 
However, Strangman and Kattamis [58] 
performed a similar experiment (at 1-g) 
and their data is provided for com- 
parison on the right side of Figure 14. 
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2,3.2 (Continued) 

Their material was of a different starting composition from the MSFC 
samples, complicating comparison. Qualitatively, agreement is good. 

The most direct omparison between theory and the MSFC experiments is 
shovm in Figure 15. Here atomic fraction is plotted versus gravity 
for a aT of 10°K at a position of z = -0.75. Although three data points 
are not really sufficient to define a curve, it is clear that the 
position and slope of any curve through the experimental points would be 
quite different from those of the theoretical curve. This difference is 
attributed to the roughly linear increase in segregation with undercooling 
exhibited by the original data (Table IV and [57]), and not accounted 
for by the AT dependence of y in equation 103. 


The hypothesis that observed Pb dendrite segregation is the result of 
nuclei transport in the melt during a transport time, t^, which is a 
function of undercooling and during which time the entire melt is under- 
cooled. Was tested by estimating the time required for a nucleus moving 
with velocity 

u = gyD [1 - (104) 

fpb 

to travel a distance of 0.75 cm. The relation of t^ to the rest of the 
solidification cycle is shown in Figure 16. For Pb nuclei of mass on 
the order of 10 gm, 


n - iiiJi 
" 2kT 


2.4 X 10"^ sec 


Thus 


- »75 cm 
t " U 


3 X 10^ 

g[i-pT^7^T 


(105) 


(106) 
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2.3.2 (Continued) 



For nominal values of g and D, the smallest values, of t. will be for the 

^ 4 

case of p{z)/ppj^ « 1. Even then, t^ is on the order of 10 sec. Since 
t is known experimentally to be about 10 sec. to 60 sec. depending on 

S 

g level and undercooling, t^ must be on the order of 10 sec. This 
discrepancy invalidates the original hypothesis, if the assumed magnitude 
of 10“^ cm^/ sec for the "transport coefficient", D, is correct for 
nuclei this large. 

Thus the niechanisnv by which undercooling produces the non-gravity- 
related segregation of Pb dendrites observed in the HSFC centrifuge 
experiments remains a mystery. In fact, the observed increase in 
segregation with increase in undercooling seen in Table IV, is contrary 
to all explanations of melt dynamics considered in this study. More 
effort will be required to explain this anomalous effect. 
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2,3.3 Morphology and Grain Size of Fe-.25 Ni Alloy 


As mentioned in Section 2,2 J ^ Kattamis an^i Flemings [56] performed the 
experimental work which established the efiect of undercooling on 
dendrite morphology and grain size. They v/orked with several alloysi but 
one which they studied most extensively was Fe-.25Ni. As so happens, 
the elements Fe and Ni are also two of the few elements, solid at STP, 
for which all the melt properties required by the mathematical model of 
Section 2.2.7 are either knovw or can be deduced. Hence this was the 
material chosen for comparison between theory and experiment. 

The first step was to evaluate all the constants required by the model 
equations. Using values of the element properties given in Table VII 
for iron and nickel, values of the constants v/ere computed from their 
defining equations. Some of the constants pertain only to the solvent 
or solute, and evaluation is straightforward. In cases v/here the 
constant pertains to the melt as a whole, the melt properties of interest, 
y, v/ere estimated from the expression 


y = ,75ypg -i- .25y„. 


007) 


where and yj^^ refer to the property for pure iron or nickel. 


Although this expression probably does not accurately represent many 
alloy melt properties, it does represent some, and it is the only means 
at present of obtaining those properties not specifically measured. The 
form of the model equations used v/as as follov/s: 


Nucleus Volume: 

Dendrite Growth Volume: 
Growth Rate: 


= I ^ f(o) = Y-, 


f(e) 

at^ 


Q o 9 ^ * H 

r*- sin 0 6t^ = V25ii> 9— 


U_ = 


vDS ^ 
c 

47)0 kT 




Nuclei Creation: 


N = I V 6t^ 


(108) 

(109) 

(no) 

( 111 ) 
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Nucleatlon Rate: 

Nucleus Radius: 

Solute Diffusion 
Constraint: 

Heat Flow Constraint; 

Grain Dimension or 
Dendrite Spacing: 


I = exp[75.985 - 

T M 


_ 2c; 

S AT 
c 






6t, 

h 

Ji = 


_ .V^2/3 /T 

“ Zaj 

(V)l/3 = 


(112) 

013) 

OH) 

015) 

(116) 


with the constants being 
32 ro \3 




IT 


(“)'" = 1 .603 X 10“^^ cm^ 
^c 

3 o./ 


s = 4 2446 X 10"^*^ 

Y2 2k A iu 


j\^ = = 7^469 X 10^ 

3 k 


M “ 4il7Pf,ia^ = 1.85019 X 10'® ^ °K 


The equations of the model (108-116) may now be solved simultaneously 
using the above constants and assumed values for the independent 
variables, e and AT. Table V was constructed to provide input values 
of f(o) and sin 6, and a-Vfang 700A computer program was written to 
facilitate calculation of several solid quantities (such as morphology 
parameter, Vg/v^^; dendrite arm spacing, 1; total dendrite volume, 

Vg + and grain size, ^) as a function of a wide range of wetting angle 
and undercooling. Table VI presents the results of such a calculation for 
the critical region for the Fe~.25N^. alloy of e=81® and 148°K < aT < 200“K, 
Kattamis and Flemings had reported a sharp transition from cylindrical 
to spherical morphology at aT =170°, hence the interest in this region. 
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TABLE V FUNCTIONS OF NUCLEANT WETTING ANGLE 


1 

0 

f(0) 

s1n^e 

0 

f(0) 

sin^o 

16 

-oonii 

.076 

56 

.12432 

.68730 

17 

.001411 

.0855 

57 

.13191 

.70337 

18 

.001767 

.0955 

58 

.13976 

.71919 

19 

.002186 

.106 

59 

.14788 

.73474 

20 

.002673 

• 

.117 

60 

.15625 

.75 

21 

.003235 

.1284 

61 

.16488 

.76496 

22 

.00388 

.14033 

62 

.17376 

.77960 

23 

.004614 

.15267 

63 

.1829 

.79389 

24 

.005444 

.16543 

64 

.19228 

.80783 

25 

.006378 

.17861 

65 

.20191 

.82139 

26 

.007423 

.19217 

66 

.21177 

.83456 

27 

.008586 

.20611 

67 

.22186 

.8473 

28 

.009875 

.2204 

68 

.23219 

.8597 

29 

.011297 

.23504 

69 

.24273 

.8716 

30 

.012861 

.25 

70 

.25349 

.883 

31 

.014572 

.26526 

71 

.26445 

,894 

32 

.01644 

.28081 

72 

.27561 

,9045 

33 

.01847 

.2966 

73 

.28697 

.9145 

34 

.02067 

.3127 

74 

.29851 

.924 

35 

.02305 

.32899 

75 

.31022 

.933 

36 

.02561 

.3455 

76 

.3221 

.9415 

37 

,02837 

.36218 

77 

.33413 

.9494 

38 

.031323 

.379 

78 

.34631 

.9568 

39 

.034481 

.396 

79 

.35863 

.9636 

40 

.03785 

.4132 

80 

.37107 

.9698 

41 

.041436 

.43041 

81 

.38363 

.9755 

42 

.04524 

.4477 

82 

.39629 

.9806 

43 

.04928 

.4651 

83 

.40905 

.9851 

44 

.05355 

.4826 

84 

.42189 

.9891 

45 

.058058 

.50 

85 

,4348 

.9924 

46 

.0628 

.51745 

86 

.44777 

.9951 

47 

.067804 

.53488 

87 

.46078 

.99726 

48 

.07305 

.5523 

88 

.47384 

.99878 

49 

.07855 

.56959 

89 

.48691 

.9997 

50 

,0843 

.5868 

90 

.5 

■ 1 1 

51 

.09032 

.60396 




52 

.0966 

.62096 




53 

.10313 

.6373 




54 

.1099 

.6545 




55 

.11699 

.67101 
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TABLE VI SOLIDIFICATION MODEL GRAIN SIZE/MORPHOLOGY 
PREDICTION 


Material; Fe - .25 Ni alloy Vfetting Angle; e = 81° 


AT 

°K 

V /v 

g n 

1 

cm 

^dendrite 

3 

V/N 

3 

ctir 

148.0 

26.1 

1.78 

5.1 X 10"^^ 

4.2 X 10"^^ 

5.68 

149.0 

21.7 

1.45 

. 

3.04 

150.0 

18.0 

1.18 

3.5 X 10"]q 
1.9 X 10 

1.65 

153.4 

9.95 

0.61 

.22 

156.0 

6.5 

.376 

1.2 X 10"t^ 
0.9 X 10"'-^ 

5.3 X 10"p 
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1.9 X 10“^ 

1S8.G 
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2.3.3 (Continued) 

In Figure 17, the theoretical prediction for morphology parameter is 

plotted beside a schematic representation of Kattamis and Flemings' 

findings. The general shape of the experimental curve is assumed to 

be similar to the theoretical, with a discontinuous change of v /v >1 

g n 

to <1 at AT = 170°K. To be consistent with Kattamis and Flemings' 

results, the experimental curve must also remain less than v /v = 10 in 

g n 

the interval 30°K < AT < 170“K, as is indicated by the onset of asymtotic 
behavior at about AT = 154°K in Figure 17. In Figure 18, the theoretical 
and experimental curves for dendrite spacing {or grain diameter) are 
compared. Although the shapes of the curves are similar, there is a 
drastic discrepancy betv/een the undercooling scales of the two graphs. 

The discrepancies found between the theoretically predicted morphological 
behavior of Fe-.25N1 alloy during solidification and that which is 
observed experimentally can be attributed to the simplifying assumptions 
made to derive the solidification modal in its present form, and to 
the almost complete lack of data for the properties of alloy melts. 

The facts that does show a transition from greater than one 

to less than one at 170°K, albeit a continuous transition, and that the 
theoretical and experimental curves for dendrite spacing are qualitatively 
similar, provide grounds for hope that an accurate predictive model for 
alloy solidification behavior (as a function of aT and e) could be 
developed with sufficient effort. The failure of the present model to 
show the proper asymtotic and discontinuous behavior of the morphology 
parameter, and its failure to match the experimental values of dendrite 
spacing (except in the region near aT = 180°K) provide an indication 
of how much more effort is needed. The usefulness of a successful and 
accurate model for predicting solid properties vn 11 be discussed in the 
next Section. 
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FIGURE 17 MORPHOLOGY PARAMETER VS. UNDERCOOLIMG FOR 
Fe-.25 N1 ALLOY, o = 81° 

. ■ 2 - 66 :,,.. 
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FIGURE 18 DEUCRITE ELEMENT SPACING VS, UNDERCOOLING FOR 
Fe-.2S Ni ALLOY, e > 81° 
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2.4 CORRELATION OF THEORY WITH EXPERIMEr'lT 

2.4.1 Recommended Procedures and Model Development 

Many procedures for the correlation of theory v/ith experiment have been 
illustrated, if not explicitly stated, in Sections 2.2 and 2,3. The 
basic problem is to express theory in a manner which predicts the value 
of parameters which are experimentally measured as a function of para- 
meters which are experimentally controlled. This usually results in 
the construction of a mathematical model, which, if the theory is 
sufficiently general, can predict the outcome of a given experiment - 
values of several parameters - as a function of those constant and 
variable parameters vfnich characterize the experimental conditions. 

A good example of this procedure is the definition of the dendrite 
morphology parameter in Section 2.2.7 and its calculation in Section 
2.3.3. Actual comparison with experimental data is often best 
accomplished with graphical techniques, though other methods such as 
tables, diagrams or theoretical versus empirical equations may be more 
informative in specific instances. 

Specific procedures will depend on the type of experiment to be per- 
formed and the nature of the theory applicable to the experiment. 

These will be defined by the overall purpose of the work. In the 
discipline of Solidification in Space, the purpose of theory and 
experiment should be to build a foundation for the creation of new 
and/or superior materials in the v;eightless environment of space, 
which will meet the requirements of the material user for certain 
specified properties possessed by the material. A methodology for 
achieving this purpose is illustrated in Figure 19, which could .be 
considered a logic diagram for a comprehensive material property 
prediction model, showing how user requirements can be matched to 
material properties by varying the input parameters which consist of 
the characteristics of the melt and its external environment - the 
solidification process controls. Although this diagram indicates 
a trial -and-error approach to selecting the pv’oper melt components 
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2.4.1 (Conti luiecl) ‘ 

and environsnent, it is orders-of “magnitude cheaper to do this analytically 
than experimentally. There also exists the possibility that as the model 
matures, methods will be devised to compute required melt components and 
environment directly from the user requirements without the need for 
trial-and-error. 

In addition to serving as a basic materials design tool, such a 
comprehensive property prediction model v/ould be invaluable to the Space 
Processing Program as an aid in selecting those materials for flight 
experiments which are most likely to produce unique or interesting properties 
in the space-processed material. Since more than half of the 92 natural 
elements are used as primary metals or as alloying elements, the number 
of possible new materials (binary, ternary, quarternary and quinary) 
is greater than 100 million. Therefore a detailed and rapid method 
of selecting materials and concentration combinations of potential 
interest is needed. ' 

The comprehensive property prediction model represented by Figure 19 
is composed of several elements, one of which is the solidification model 
described in Section 2.2.7. Development of each element of the com- 
prehensive model will generate a number of vc”if1 cation experiments, in 
both one-g and zero-g, to confirm that model predictions correspond to 
nature, particularly that part of nature which is so removed from 
everyday experience: zero-g. The concept of a comprehensive material 

property prediction model is relatively new, and at present only the 
solidification model mentioned above and the melt property input for 
the Fe-Mi and Pb-Sn systems have received much attention. But from these, 
certain procedures for correlating theory with experiment can be 
ascertained. 

The primary factor in determining what materials are to be used in the 
comparison of theory v;ith experiment is knowledge of the melt properties 
of the chosen material. If such properties as thermal diffusivity, 
entropy of fusion, liquid-solid interfacial energy, diffusion coefficient 
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2.4.1 (Continued) 

and density are not known, then the input constants for the solidification 
model cannot be specified, and nothing can be calculated for the material. 
Hence, those materials for which melt property values are known should 
be prime candidates for experiment. Other factors, such as the 
benefit of a large density differential in the Pb-Sn system, must also 
be considered in the material selection process. 

After selection of materials and definition of the input variables and 
constants, the quantities to be measured experimentally and predicted 
theoretically must be determined. Then building of the mathematical 
model and design of the experiment may proceed. If the experiment has 
already been performed, as in the case of the Fe-.25 Hi morphology 
study in one-g, then the quantities to be predicted are specified by 
the experiment- If the experiment has not been designed, then theory will 
provide a guide as to what quantities are most desirable or practicable 
to measure. After performing the experiment, making measurements and 
carrying out computations with the model, the data obtained should be 
compared by plotting it on the same graph, displaying it in the same 
table, etc. In this way the accuracy of the theory in predicting or 
explaining physical reality can be determined. More specific discussion 
of comparison procedures and experiment objectives v/ill be given in the 
following Sections. It is recommended that NASA begin, as soon as 
possible, to develop a comprehensive material property prediction model, 
similar to that illustrated by Figure 19, to serve as a theoretical 
guide for Space Processing development, and later as a molecular-level 
materials design tool. 

2.4.2 Recommendations of Materials for Flight Experiments 

As mentioned above, a primary consideration in selecting a material 
system for study is the availability of data for the liquid state 
properties of the material. Even v/hen data exists for a material, it 
may be very difficult to obtain in a form useful for solidification 
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2,4.2 (Continued) 

calculations. A case in point is thermal diffusivity, ay. Usually 
thermal diffusivity values are not reported in the literature for 
liquids, and ay must be calculated from its defining equation 

K 


(117) 


where k: is thermal conductivity and Cp is specific heat. The melt 

values of k, p and c^ must now be obtained. Values of p and c„ for 

P P 

liquid elements are fairly common, but fc must often be calculated from 

electrical conductivity, a^, data via the W1edman-Franz law 


3e^ 


(118) 


v\'here e is the electron charge. 

In addition to the almost complete lack of experimentally determined 
values for the properties of multi -component melts, there is a dearth 
of techniques for calculating melt property values for alloys. When 
all else fails, the simple relation 


^melt 


(119) 


v;here A is the solute atomic fraction and y^. is the melt value of property 
y for component i, may be used to obtain an order^-of-magnitude estimate 
for the Pi'operty However, with certain exceptions [59], there 

is no rigorous basis for using equation 119 to obtain alloy melt 
properties. Other criteria for selecting a material system for theoretical 
and experimental study include: 


6 Availabiiity of one-g experiment data bank, 

G The practical or scientific importance of the material system, 

9 The properties of. the material system enhance a particular measurement 
to be made, e.g., different component densities enhance segregation 
measurement, 
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2.4.2 (Continued) 

D Simplicity of the material system with reference to experiment 
design and theoretical analysis. 

Taking these criteria into account ^ the two material systems recommended 
for further study, particularly zero-g experimentation, are iron-nickel 
and lead-tin. Melt property data for the elemental components of these 
systems is given in Table VII, along with solute diffusion coefficients. 
Alloy melting temperatures are given by phase diagrams, the atomic 
weight, density and specific volume can be found by equation 119, and 
the remaining alloy properties are unknown. 

2.4.3 Recommended Experiment Objectives 

The primary purpose of the study described by this report has been to 
define theoretically how undercooling, gravity and nucleants (impurities, 
high temperature phases or container surfaces) affect the solidification 
of materials and the resultant solid properties. Hence it is recommended 
that zero-g experiments (with one-g controls) be designed v/hich will 
determine the effects of these environmental variables on solidification 
and solid properties. Several specific experiment types can be identified: 

Segregation Experiments 

These experiments vvouTd essentially repeat the MSEC centrifuge experiments 
vnth Pb-Sn eutectic melts, but be performed in zero-g to isolate the 
unexplained undercooling effect on solute segregation from the gravity 
effects. Volume fraction (or atomic fraction) of lead and temperature 
(undarcpol ing) would be the primary variables measured. 

Morphology Experiments 

Morphology experiments would be designed to measure nucleati on rate, 
growth rate, grain size and dendrite morphology as functions of under- 
cooli ng and nucTeant wetti ng angl e (or other parameter appropriate to 
characterization of the nucleant). The undercooling and rate measurements 
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TABLE VII ■ LIQUID STATE/SOLIOIFICATION PARAMETERS FOR 
SELECTED MATERIALS: 


^T 

m 

M 

gm 

gm atom 

p 

9”/cm3 

TJ 

3 

cm 

0 ^ 
ergs 
cm^ 

ergs 

K 

watt 

cm 

Cp 

joule 

gm°K 

■ 

cm^ 

sec 

gm atom 


Fe 

1809 

55,85 

6.9 

8.09 

204 

12X10^ 

0.39 ■ 

0.7866 

0.072 

Ni 

1725 

i 

58,71 

8,0 

■ 7.34 

255 

15.5X10^ 

0.469 

0.6556 

0,089 

Pb 

600,6 

207.21 

i 

10.3 

i 20.1 

33.3 

4.34X10^ 

0.165 

1 

■ 1 

! 

0.1417 

0.106 

Sn 

505 

118.7 

6.98 

17.0 

59 

8.53X10^ 

0.369 

0.2590 

0.216 

REFERENCES 

PERIODIC 

PERIODIC 

INT. CRIT. 

M/ 

D. Tu: 

RNBULL 

INT. CRIT. 

HOLTGREN, 

K 


TABLE 

TABLE 

TABLES 

/p 



TABLES 

et . al . 



2 

DIFFUSION DATA: ^^Ni—iFe = 4,67X10'^ /sec 

2 

Dpb_^Sn = 3.69X10"^ /sec 
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v/Quld be made during the zero-g flight period of the experiment^ and 
morphological measurements could be made later after sectioning the 
sample. In making morphological measurements j a method should be 
devised for measuring the average size or volume of the dendrites. 

From the dendrite volum-.’, v^, the morphology parameter may be obtained 
directly by 

V''n = V^'n - V 

vdiere the nucleus volume » V|^, is obtained from the undercooling and 
watting angle by equation 95. 

Undercooling and Nucleation Experiments 

The purpose of these experiments would be to determine the effects of 
gravity and nucleants, including container surfaces, on melt undercooling. 
Of particular interest are convection effects on undercooling (temperature) 
distribution in the melt and the influence this distribution has on . 
nucleus formation and growth-especially on interface stability. A goal 
of containerless solidification experiments in zero-g, is to establish 
the importance of container surfaces as nucleants which limit the 
achievable, undercooling, and thus, control over solid properties. The 
control experiment in this case would be another zero-g sample v/hich 
2 ^ in a container. In all cases the measurement required vn 11 be 
undercooling (temperature) as a function of position in the melt. 

Additional measurements of convective heat transfer and g level v;ill be 
desirable in the gravity-dependent experiments. 

M ud eant Characteri zati on Experiments 

These experiments will attempt to specify the most convenient parameter 
for characterizing the influence v/hich nueleant interfacial energies have 
on the nucTeation and growth processes. Paranieters used in this report 
include nucleus-melt interfacial energy, o (equation 16) and 

Ao (equations 21). If, for example, it is found that most nuclei 
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2.4.3 (Continued) 

which form are of the spherical cap type, a further purpose of these 
experiments would be to correlate 0 values v/ith particular impurity 
(nucleant) - melt combinations and to hopefully provide insight into 
the mechanism \^hich causes a nucleus to form more readily on one 
type of nucleant than on another. If this insight is obtained, it 
should provide some understanding of denucleation as v/ell . 

Melt Property Measurements 

As mentioned previously, the main problem encountered in the calculation 
of theoretical predictions of solidification phenomena (disregarding the 
accuracy of the theory itself), is the dearth of values for the properties 
of the melt V'/hich usually appear as constants in the equations of the 
theory. Since most of these -properties are defined as fundamental constants 
for the material or are derived from intermblecular structure and energy 
in a formal manner which can not yet be reduced to numerical calculation, 
experimental measurement of these properties is the most direct method for 
obtaining them. This is usually quite difficult for metals and alloys 
due to the high temperatures required by the liquid state and the sometimes 
corrosive nature of the metal melt being studied. But if the ability to 
predict the results of solidification processes and to enter an era of true 
design of material properties at the molecular level is to be attained, 
values of melt properties such as thermal diffusivity, interfacial energy 
(or a corresponding parameter), diffusion coefficient, specific volume 
(density) , i sothermaT compressibility, thermal expansion coeffi cient 
and entropy of fusion must be obtained. 
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